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ABSTRACT 

The geometry of the cosmic web drives in part the spin acquisition of galaxies. This can 
be explained in a Lagrangian framework, by identifying the specific long-wavelength 
correlations within the primordial Gaussian random field which are relevant to spin 
acquisition. Tidal Torque Theory is revisited in the context of such anisotropic en¬ 
vironments, biased by the presence of a filament within a wall. The point process of 
filament-type saddles represents it most efficiently. The constrained misalignment be¬ 
tween the tidal and the inertia tensors in the vicinity of filament-type saddles simply 
explains the distribution of spin directions. This misalignment implies in particular 
an azimuthal orientation for the spins of more massive galaxies and a spin alignment 
with the filament for less massive galaxies. This prediction is found to be in qualitative 
agreement with measurements in Gaussian random fields and N-body simulations. It 
relates the transition mass to the geometry of the saddle, and accordingly predicts its 
measured scaling with the mass of non-linearity. Implications for galaxy formation and 
weak lensing are briefly discussed, as is the dual theory of spin alignments in walls. 

Key words: cosmology: theory — galaxies: evolution — galaxies: formation — galax¬ 
ies: kinematics and dynamics — large-scale structure of Universe — 


1 INTRODUCTION 

Modern simulations based on a well-established paradigm of 
cosmological structure formation predict a significant con¬ 
nection between the geometry and dynamics of the large- 
scale structure on the one hand, and the evolution of the 
physical properties of forming galaxies on the other. Key 
questions formulated decades ago are nevertheless not fully 
answered. What are the main processes which determine the 
morphology of galaxies? What is the role played by angular 
momentum in shaping them? 

Pichon et al. (2011) have suggested that the large-scale 
coherence of the inflow, inherited from the low-density cos¬ 
mic web, explains why cold flows are so efficient at producing 
thin high-redshift discs from the inside out (see also Stew¬ 
art et al. 2013; Laigle et al. 2015; Prieto et al. 2014). On the 
scale of a given gravitational patch, gas is expelled from ad¬ 
jacent voids, towards sheets and filaments forming at their 
boundaries. Within these sheet s/filaments, the gas shocks 
and radiatively loses its energy before streaming towards 
the nodal points of the cosmic network. In the process, it 
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advects angular momentum, hereby seemingly driving the 
morphology of galaxies (bulge or disc). The evolution of the 
Hubble sequence in such a scenario is therefore at least in 
part initially driven by the geometry of the cosmic web. As 
a consequence, the distribution of the properties of galaxies 
measured relative to their cosmic web environment should 
reflect such a process. In particular, the spin distribution of 
galaxies should display a preferred mass-dependent orienta¬ 
tion relative to the cosmic web. 

Both numerical (e.g. Aragon-Calvo et al. 2007; Hahn 
et al. 2007; Sousbie et al. 2008; Paz et al. 2008; Zhang et al. 
2009; Codis et al. 2012; Libeskind et al. 2013; Aragon-Calvo 
2013; Dubois et al. 2014), and observational evidence (e.g. 
Tempel et al. 2013) have recently supported this scenario. 
In parallel, much analytical (e.g. Catelan et al. 2001; Hirata 
& Seljak 2004), numerical (e.g. Heavens et al. 2000; Croft & 
Metzler 2000; Schneider & Bridle 2010; Schneider et al. 2012; 
Joachimi et al. 2013b; Codis et al. 2015; Tenneti et al. 2015) 
and observational (e.g. Brown et al. 2002; Lee & Pen 2002; 
Bernstein & Norberg 2002; Heymans et al. 2004; Hirata et al. 
2004; Hirata Sz Seljak 2004; Mandelbaum et al. 2006; Hirata 
et al. 2007; Mandelbaum et al. 2011; Joachimi et al. 2011, 
2013a) efforts have been invested to control the level of in- 
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trinsic alignments of galaxies as a potential source of sys¬ 
tematic errors in weak gravitational leasing measurements. 
Such alignments are believed to be a worrisome source of 
systematics of the future generation of lensing surveys like 
Euclid or LSST. It is therefore of interest to understand from 
first principles why such intrinsic alignments arise, so as to 
possibly temper their effects. 

Hence we should try and refine a theoretical frame¬ 
work to study the dynamical influence of filaments on 
galactic scales, via an extension of the peak theory to the 
truly three-dimensional anisotropic geometry of the circum- 
galactic medium, and amend the standard galaxy formation 
model to account for this anisotropy. Toward this end, we 
will develop here a filament version of an anisotropic “peak- 
background-split” formalism, i.e. make use of the fact that 
walls and filaments are the interference patterns of primor¬ 
dial fluctuations on large scales, and induce a corresponding 
anisotropic boost in over-density. Indeed, filaments feeding 
galaxies with cold gas are themselves embedded in larger 
scale walls imprinting their global geometry (Danovich et al. 
2012; Dubois et al. 2012). On top of these modes, construc¬ 
tive interferences of high frequency modes produce peaks 
which thus get a boost in density that allows them to pass 
the critical threshold necessary to decouple from the overall 
expansion of the Universe, as envisioned in the spherical col¬ 
lapse model (Gunn & Gott 1972). This well-known biased 
clustering effect has been invoked to justify the clustering of 
galaxies around the nodes of the cosmic web (White et al. 
1988). It also explains why galaxies form in filaments: in 
walls alone, the actual density boost is typically not suffi¬ 
ciently large to trigger galaxy formation. The main nodes 
of the cosmic web are where galaxies migrate, not where 
they form. They thus inherit the anisotropy of their birth 
place as spin orientation. During migration, they may col¬ 
lide with other galaxies/haloes and erase part of their birth 
heritage when converting orbital momentum into spin via 
merger (e.g. Godis et al. 2012). Tidal torque theory should 
therefore be re-visited to account for the anisotropy of this 
filamentary environment on various scales in order to model 
primordial and secondary spin acquisition. 

In this paper, we will quantify and model the intrinsi¬ 
cally 3D geometry of galactic spins while accounting for the 
geometry of saddle points of the density field. Indeed, sad¬ 
dle points define an anisotropic point process which accounts 
for the presence of filaments embedded in walls (Pogosyan 
et al. 1998), two critical ingredient in shaping the spins of 
galaxies. Taking them into account will in particular al¬ 
low us to predict the biased geometry of the tidal field in 
the vicinity of saddle points. This can be formalized using 
the two-point joint probability of the gravitational poten¬ 
tial field and its first to fourth derivatives and imposing a 
saddle point constraint. For Gaussian (or quasi-Gaussian) 
helds these two-point functions are within reach from first 
principle (Bardeen et al. 1986). A proper account of the 
anisotropy of the environment in this context will allow us 
to demonstrate why the spin of the forming galaxies field are 
first aligned with the filament’s direction. We will also show 
that massive galaxies will have their spin preferentially along 
the azimuthal direction. While relying on a straightforward 
extension of Press Schechter’s theory, we will predict the 
corresponding transition mass’ scaling with the (redshift de¬ 
pendent) mass of non-linearity, while relying on the so-called 



Figure 1. Spin acquisition by tidal torquing. At linear order, the 
misalignment between the inertia tensor of the proto-object and 
the surrounding tidal tensor induces an inhomogeneous ZePdovich 
boost which corresponds to the acquisition of a net intrinsic an¬ 
gular momentum in Eulerian space. 

cloud-in-cloud problem applied to the filament-background 
split. 

The paper is organised as follows. Section 2 qualita¬ 
tively presents the basis of the physical process at work in 
aligning the spin of dark haloes relative to the cosmic web. 
Section 3 then presents the expected Lagrangian spin distri¬ 
bution near filaments, assuming cylindrical symmetry, and 
explains the observed mass transition while carrying a multi¬ 
scale analysis of the fate of collapsing haloes in the vicinity 
of 2D saddle points. Section 4 revisits this distribution in 
three dimensions for realistic typical 3D saddle points. Sec¬ 
tion 5 investigates the predictions of the theory using Gaus¬ 
sian random field and N-body simulations, while we finally 
conclude in Section 6. Appendix A discusses possible limita¬ 
tions and extensions of this work. Appendix B presents the 
dual theory for spin alignment near wall saddles. Finally, 
Appendix G gathers some technical complements. 


2 TIDAL TORQUING NEAR A SADDLE 

Before presenting analytical estimates for the expected spins 
near filament in two and three dimensions and their tran¬ 
sition mass, let us discuss qualitatively what underpins the 
corresponding theory. 

2.1 Spin acquisition by tidal torquing 

In the standard paradigm of galaxy formation, protogalax¬ 
ies acquire their spin^ by tidal torquing coming from the 
surrounding matter distribution (Hoyle 1949; Peebles 1969; 
Doroshkevich 1970; White 1984; Gatelan & Theuns 1996; 
Grittenden et al. 2001). At linear order, this spin is acquired 
gradually until the time of maximal extension (before col¬ 
lapse) and is proportional to the misalignment between the 
inertia tensor of the protogalaxy and the surrounding tidal 
tensor (see Schaefer 2009, for a review) 

Li = ^ (t)D-^(t)eijkIjiTik , (1) 

j,k,l 

where a{t) is the scale factor, D+ the growth factor, Tij the 
tidal tensor (detraced Hessian of the gravitational poten¬ 
tial), lij the protogalactic inertia tensor (only its traceless 

^ Note that in this paper we will call interchangeably “spin” or 
“angular momentum” the intrinsic angular momentum of (proto-) 
haloes. 
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part, lij contributes to the spin). As this work focuses on 
the spin direction, the factor a^(t)D+(t) will henceforth be 
dropped for brievity. This process of spin acquisition by tidal 
torquing is illustrated on Figure 1. 

In the Lagrangian picture, Uj is the moment of inertia 
of a uniform mass distribution within the Lagrangian image 
of the halo, while Tij is the tidal tensor averaged within the 
same image. Thus, to rigorously determine the spin of a halo, 
one must know the area from which matter is assembled, 
beyond the spherical approximation. While this can be de¬ 
termined in numerical experiments, theoretically we do not 
have the knowledge of the exact boundary of a protohalo. As 
such, one inevitably has to introduce an approximate proxy 
for the moment of inertia (and an approximation for how 
the tidal field is averaged over that region). 

The most natural approach is to consider that proto¬ 
haloes form around an elliptical peak in the initial density 
and approximate its Lagrangian boundary with the ellipti¬ 
cal surface where the over-density drops to zero. This leads 
to the following approximation for the traceless part of the 
inertia tensor (e.g. Schafer & Merkel 2012, see also equa¬ 
tions (A2)-(A4)) 

- 2 2 M ~ 

lij = 

where is the traceless part of the inverse Hessian of the 
density field, Hij = didjS, v is the overdensity at the peak, 
and M is the mass of the protohalo. In the second form we 
explicitly presented the inverse Hessian via the (detraced) 
matrix of the Hessian minors, Hij. While Hij is a simple 
polynomial in second derivatives of the density, M/detH is 
not, which is the source of most technical difficulties when 
statistical studies of the spin are attempted. 

Let us point at the following considerations to bypass 
these difficulties. First, the supplementary condition for the 
halo to be at a peak of the density yields an extra det H fac¬ 
tor in all statistical measures (see e.g. Bardeen et al. 1986). 
This factor exactly cancels the determinant in the denomi¬ 
nator. Secondly, all quantities in equation (2) are computed 
after the density field is smoothed at a particular scale Rh 
which sets the corresponding mass scale. Therefore, it is 
more appropriate to apply equation (2) to haloes at fixed 
mass M, determined by that smoothing. Hence we could ar¬ 
gue for the proxy Rj oc vHij for the moment of inertia for 
haloes of a given hxed mass, where the change in mass is re¬ 
flected in the corresponding change in the smoothing scale. 
In two dimensions, we show in Appendix A that this multi¬ 
scale approximation gives qualitatively the same statistical 
results as just using Hij as a proxy. While this approxima¬ 
tion is relatively simple, since we are only concerned with 
the direction of the spin, we will now go one step further 
and use throughout this paper the Hessian as a proxy for 
the inertia tensor, even in three dimensions. Indeed, Aj, Hij 
and Hij share the same eigen-directions (Catelan & Theuns 
1996; Schafer & Merkel 2012), so we define the spin for the 
rest of the paper as 

Si = ^ ^ CijkHjlTlk, . (3) 

j,k,l 

The vector field Si is then quadratic in the successive deriva¬ 
tives of the potential: its (possibly constrained) expectation 
can therefore be computed for Gaussian random fields. This 



Figure 2. On top of the density contours (from dark blue to 
dark red), the three (red triangle) maxima (resp. the three (green 
points) minima) are connected by the crest lines (in solid gold, 
resp. the through-lines in dashed green) which intersect through 
saddles points (purple squares). The blue arrows represent stream 
lines of the gradient flow. Throughout this paper, we will assume 
that the geometry of the regions of intermediate densities are set 
by the shape of the (purple) saddles. 


approximation is further discussed in Appendix A. Note that 
equation (3) improves upon simple parametrisations of the 
mean misalignment between inertia and tidal tensors (see, 
e.g. Lee & Pen 2000; Crittenden et al. 2001) by ah initio 
explicitly taking into account the correlations between both 
tensors. 


2.2 Geometry of the cosmic web 

Galaxies are not forming everywhere but preferentially in 
filaments and nodes which define the so-called cosmic web 
(Klypin & Shandarin 1993; Bond et al. 1996). The origin of 
these structures lies in the asymmetries of the initial Gaus¬ 
sian random field describing the primordial universe, ampli¬ 
fied by gravitational collapse (Zel’dovich 1970). The pres¬ 
ence of such large-scale structure (walls, filaments, nodes) 
induces local preferred directions for both the tidal tensor 
and the inertia tensor of forming objects which will even¬ 
tually turn into preferred alignments of the spin w.r.t the 
cosmic web. It is therefore of interest to understand what is 
the expected spin direction predicted by equation (3) given 
the presence of a typical hlament nearby. As a filament is 
typically the field line that joins two maxima of the den¬ 
sity field through a filament-type saddle point (where the 
gradient is null and the density Hessian has two negative 
eigenvalues), we choose to study in this paper the expected 
spin direction of proto-objects in the vicinity of a filament- 
type saddle point with a given geometry (which imposes the 
direction of the filament and the wall). Figure 2 illustrates 
the geometry of filaments near peaks and saddles in a 2D 
Gaussian field. 
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2.3 Constrained tidal torque theory in a nutshell 

2.3.1 spin alignments and flips 

It has been shown in simulations (Bailin & Steinmetz 2005; 
Aragon-Calvo et ah 2007; Paz et ah 2008; Zhang et ah 2009; 
Codis et ah 2012; Libeskind et ah 2013; Forero-Romero et ah 
2014, among others) that the spin of dark haloes is corre¬ 
lated to the direction of the filaments of the cosmic web in 
a mass-dependent way. The alignment between the spin and 
the closest filament increases with mass until a mass of max¬ 
imum alignment (Laigle et ah 2015) that we call here critical 
mass. As mass increases, the direction of the spin becomes 
less aligned with the filament before becoming perpendicu¬ 
lar to it (Codis et al. 2012). This transition - from aligned 
to perpendicular - occurs at a mass that we call here the 
transition mass. 

This paper will claim that the critical mass is directly 
related to the size of the quadrant of coherent angular mo¬ 
mentum imposed by the tides of the saddle point (which 
are effectively the Lagrangian counter parts of the quad¬ 
rant of vorticity found in Laigle et al. (2015)). This mass 
can be captured using a cylindrical model that would cor¬ 
respond to the plane perpendicular to the filament at the 
saddle point (which amounts to assuming an inhnitely long 
filament). This 2D toy model (see Section 3 below) shows 
that near a 2D peak (i.e near an infinitely long 3D filament), 
the quadrupolar structure seen in simulation naturally arises 
in a Lagrangian framework. We investigate the size of that 
quadrants and shows that it qualitatively predicts the right 
critical mass. 

The second stage of accretion, that flips the spin of more 
massive haloes from aligned to perpendicular to the fila¬ 
ments, requires a 3D analysis (see Section 4). It is shown 
that indeed small-haloes that form close to the saddle 
point, acquire spin along the filaments while more massive 
haloes that form further from the saddle (i.e closer to the 
peaks/nodes) acquire a spin perpendicular to the filaments 
(while accreting smaller haloes). The transition mass will be 
predicted as a function of redshift and shown to agree with 
measurements in simulations. 

2.3.2 The premiees of anisotropie tidal torque theory 

Let us present here an outline of the extension of TTT within 
the context of a peak (or saddle) background split. Given 
the an-isotropically triaxial saddle constraint, we will argue 
that the misalignment between the tidal tensor and the hes¬ 
sian of the density field simply explains the transverse and 
longitudinal antisymmetric geometry of angular momentum 
distribution in their vicinity. It arises because the two ten¬ 
sors probe different scales: given their relative correlation 
lengths, the hessian probes more directly its closest neigh¬ 
bourhood, while the tidal field, somewhat larger scales. 

Within the plane of the saddle point perpendicular to 
the filament axis (the midplane hereafter), the dominant 
wall (corresponding to the longer axis of the cross section 
of the saddle point) will re-orient more the Hessian than 
the tidal tensor, which also feels the denser, but typically 
further away saddle point, see Figure 3, top panels. This net 
misalignment will induce a spin perpendicular to that plane 
i.e along the filament. This effect will produce a quadrupolar, 
antisymmetric distribution of the longitudinal component 


of the angular momentum which will be strongest at some 
four points, not far off axis. Beyond a couple of correlation 
lengths away from those four points, the effect of the tidal 
field induced by the saddle point will subside, as both tensors 
become more spherical. 

Conversely, in planes containing the filament, e.g. con¬ 
taining the main wall, a similar process will misalign both 
tensors. This time, the two anisotropic features differentially 
pulling the tensors are the filament on the one hand, and the 
density gradient towards the peak on the other. The net ef¬ 
fect of the corresponding misalignment will be to also spin 
up haloes perpendicular to that plane, along the azimuthal 
direction, see Figure 3, bottom panels. By symmetry, the 
anti-clockwise tidal spin will be generated on the other side 
of the saddle point. 

Hence the geometry of angular momentum near 
filament-saddle points is the following: it is aligned with 
the filament in the median plane (within four antisymmet¬ 
ric quadrants), and (anti-)aligned with the azimuthal direc¬ 
tion away from that plane. The stronger the triaxiality the 
stronger the amplitude. Conversely, if the saddle point be¬ 
comes degenerate in one or two directions, the component 
of the angular momentum in the corresponding direction 
will vanish. For instance, a saddle point in the middle of 
a very long hlament will only display alignment with that 
filament axis, with no azimuthal component. For a typical 
triaxial configuration, two pairs of four points define the loci 
of maximal longitudinal and azimuthal spin. 

2.3.3 Geometry of spin flip 

Figure 4 gives a more quantitative account of the geometry 
of the tidal field around a given saddle point embedded in 
a given dominant wall. We consider here the angular mo¬ 
mentum distribution near a hlament-saddle point, S. It is 
assumed that the three eigenvalues of the density are such 
that the filament going through this saddle point is along the 
vertical axis and that the other two eigenvalues are differ¬ 
ent, reflecting the presence of a dominant wall, in the Sxz 
plane, in which the filament is embedded. The shape of a 
given triaxial iso-density is shown in purple, together with 
two cross sections, resp. in the Sxy mid plane and in a plane 
containing the maxima of transverse angular momentum. As 
we will demonstrate later, the spin is mostly conhned in the 
neighbourhood of the Sz axis, up to a couple of correla¬ 
tion length of the density. It would in fact vanish, should 
the saddle become isotropic. In the Sxy plane, we identify 
four quadrants corresponding to regions in which the spin 
is parallel to the filament. Within theses quadrant, the spin 
point respectively along in the the hrst and third quad¬ 
rants, and along —in the second and fourth quadrant. By 
symmetry, the spin has to vanish along Sx and Sy. 

2 . 3.4 Towards a transition mass? 

The twisted geometry of the spin near the saddle point also 
allows us to identify the Lagrangian transition mass corre¬ 
sponding to the alignment of dark matter haloes’ spin rela¬ 
tive to the direction of their neighbouring hlament. Let us 
hrst consider the Lagrangian counterpart of a low-mass halo 
and assume it lies near the median plane, it will typically 
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Figure 3. Sketch of main differential alignment between halo shapes and tidal tensor responsible for and component of momentum. 
Top: the two tensors in light and dark red, end up being misaligned as they feel differently the neighbouring wall (blue) and filament 
(purple), inducing a spin parallel to the filament (red arrow). Three projections are shown for clarity. Bottom: correspondingly, the 
differential pull from the filament (purple) and the density gradient towards the peak (blue) generates a spin (red arrow) along the 
azimuthal direction. By symmetry, the other peak(s) on the other side of the saddle point will spin up massive haloes in the opposite 
direction. 


fall into one of the quadrant corresponding to an orienta¬ 
tion of the spin parallel to the filament axis. Now consider 
a halo of larger Lagrangian extent. As long as its size is 
smaller than the typical size of a quadrant (which will be 
defined more precisely below) the alignment increases, until 
it over-extends the quadrant. As it does, two things happen 
i) it will start capturing tides from the next quadrant, which 
would anti-align it; as the Lagrangian patch radius increases 
more, it reaches a size comparable to the whole tidal region 
of influence of the saddle point. It then encompasses both 
the clock-wise and anti-clockwise azimuthal regions, and add 
up to a net momentum of null amplitude, ii) it will start 
capturing the effect of the azimuthal tide, hence inducing a 
spin flip. Depending on the ratio of the eigenvalues of the 
Hessian, the two might be concurrent or not. In parallel, as 
the radius increases, the patch collects the mean potential 
gradient which defines the Zel’dovich boost which will drive 
it away from the neighbourhood of the saddle point. The 
above description clearly accounts for the influence of only 
one saddle point. As we consider regions further away from 
that saddle, we should account for the influence of other 
critical points, as discussed in Section 5. 

We have up to now considered a patch centered near the 
midplane close to the saddle point. Indeed, typically, in the 
peak background split framework, such patches will collapse 
preferentially where the density is boosted, that is within the 
wall containing the filament, close to the filament. The rarer 
(more massive) haloes will form in turn in the denser regions, 
away from the saddle point, along the filament, while the 
more common lighter haloes will form everywhere and in 
particular near the saddle point. The former will have a 


spin perpendicular to the filament. The latter will have a 
spin parallel to the filament. The relative number of light to 
small haloes will depend on curvilinear coordinate along the 
filament because consumption is important: object above 
the transition mass have swallowed their lighter parents. At 
a given redshift, the left overs will decide what matters. This 
effect is the anisotropic version of the well-known cloud-in¬ 
cloud problem. 

2.3.5 Lagrangian dynamics of spin flip 

In order to understand how a given halo flips, let us split the 
original Lagrangian patch in two concentric shells. The inner 
shell will correspond to the Lagrangian extent of the halo 
as it initially forms, while the outer shell will correspond to 
secondary infall. The reasoning presented in Section 2.3.2 
can be applied independently to both the inner and outer 
shell, and we would typically conclude that the outer shell 
would be more likely to have its spin perpendicular to the 
filament axis. It follows that, as far as this halo is concerned, 
it will undergo a spin flip as it moves towards the core of 
the filament and away from the saddle. This process will also 
correspond to an acquired net helicity for the secondary in¬ 
fall, which will last as long as the transverse anisotropy of 
the saddle point correlates the local tidal field. In effect, 
this consistent helicity will build up the spin of the form¬ 
ing galaxy via secondary infall as it drifts, up to the point 
where mergers will re-orient the direction of its spin. Hence 
this constructive build up of disc should only last so long as 
the galaxy drifts within the high-helicity region. Note that 
the transverse motion will correspond to the halo entering 
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Figure 4. Qualitative geometry of the angular momentum dis¬ 
tribution near a elliptical saddle point S (see also Figure 12). 
The shape of a given tri-axial iso-density is shown in purple, to¬ 
gether with two cross sections, resp. in the Sxy mid plane and 
in a plane containing the maxima of transverse angular momen¬ 
tum. The velocity flow in the Sxz plane is shown in cyan. The 
locus of the Lagrangian extent of haloes is shown in orange con¬ 
centric spheres centered on M. In the Sxy plane, the four points, 
Pz correspond to the maxima of the modulus of angular momen¬ 
tum. They point respectively along in the the first and third 
quadrants (in red), and along —in the second and fourth quad¬ 
rant (in blue). Conversely, the four points, correspond to the 
maxima of the alignment of s along Only two (above the mid 
plane) are shown. As the orange sphere capture more than one 
quadrant, the 2 : component of L subside, as it encompasses more 
of the neighbourhood of P^^, and its component increases. In 
this Lagrangian framework, the spin flip as a function of mass is 
a direct consequence of the geometry of the tidal held imposed 
by the saddle point. 


the vortex rich caustic corresponding to the multi-flow re¬ 
gion near the filament, so that this Lagrangian description 
remains fully consistent with the Eulerian discussion given 
in Laigle et al. (2015). We can anticipate that the longitu¬ 
dinal motion generates azimuthal vortices as well. 

The scenario described in this section can be formal¬ 
ized at two levels. First, within the framework of constrained 
random fields, one can compute the expected geometry of 
the spin configuration near a given saddle. This will yield 
a map of the mean alignment between spin and filament in 
the vicinity of the saddle point. We will then marginalize 
over the expected distribution of such saddles, and model 
correspondingly the evolution of the expected mass of dark 
haloes around the filament. This will allow us to recover the 
numerically measured mass transition for spin flip. We may 
also test the mass-dependent alignment w.r.t. in Gaus¬ 
sian random fields and N-body simulations. For the sake 
of clarity, we will proceed in two steps: first, while assuming 
cylindrical symmetry we will compute the expected spin dis¬ 
tribution within the most likely cross section of a filament of 
infinite extend (Section 3); then we will compute this expec¬ 
tation around the most likely 3D saddle point (Section 4). 


3 SPIN ALONG INFINITE FILAMENT 

Let us first start while assuming that the filament is of in¬ 
finite extent, so that we can restrict ourselves to cylindrical 
symmetry in two dimensions. This is of interest as the angu¬ 
lar momentum is then along the filament axis by symmetry 
and its derivation in the context of TTT is much simpler. 
It captures already in part the mass transition, in as much 
as we can define the mean extension of a given quadrant of 
momentum with a given polarity. In this context, it is of in¬ 
terest to study the spin geometry in the median plane i.e in 
the vicinity of a 2D peak. This 2D spin is along the filament, 
and will be denoted in what follows. 

3.1 Shape of the spin distribution near filaments 

Under the assumption that the direction of the spin along 
the z direction is well represented by the fully anti¬ 
symmetric (Levi Civita) contraction of the tidal tensor and 
density hessian given by equation (3) (e.g. Schafer & Merkel 
2012), it becomes a quadratic function of the second and 
fourth derivatives of the potential. As such, it becomes pos¬ 
sible to compute expectations of it subject to its relative 
position to a peak with a given geometry (which would cor¬ 
respond to the cross section of the filament in the midplane). 
Note that, as mentioned in Section 2.1, standard TTT relies, 
more correctly, on the inertia tensor in place of the Hessian. 
Even though they have inverse curvature of each other, their 
set of eigen-directions are locally the same, so we expect the 
induced spin direction- which is the focus of this paper, to 
be the same, so long as the inertia tensor is well described 
by its local Taylor expansion. 

3.1.1 Constrained joint PDF near peak 

Any matrix of second derivatives fij - rescaled so that 
((A/)^) = I- can be decomposed into its trace A/, and 
its detraced components in the frame of the separation 

T = (/i1-/22)/2, /^=/i2. (4) 

Then all the correlations between two such matrices, fij and 
Qij can be decomposed irreducibly as follows. Let us call 
and the correlation functions in the frame of 
the separation (which is the first coordinate here) between 
the second derivatives of the field / and g separated by a 
distance r: 

Cf/(r) = (A/As>, 

qUw=(A/T), 

( 5 ) 

All other correlations are trivially expressed in terms of the 
above as 

(/><Ag) = 0, (/V) = 0, 

(/V) = (6) 

Here, we consider two such fields, namely the gravitational 
potential $ and the density 6 . In the following these two 
fields and their first and second derivatives are assumed to 
be rescaled by their variance ctq = 
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al = {(5 = A<I>)2), cr| = {{y5f) and cr| = ((ASf). The 
shape parameter of the density field is defined as 

-i ^ gIKp^cj^). (7) 


The rescaled potential and density will be denoted by 0 and 
X and the rescaled first and second derivatives by , Xi and 

(pij ? Xij. 

Let us gather the first and second derivatives of the 
gravitational field and the first and second derivatives of 
the density in a vector denoted by X spatially located in 
rx and Y located in ry- The Gaussian joint PDF of X and 
Y at the two given locations (rx and ry separated by a 
distance r = \rx — ry\) obeys 


P(X,Y) 


1 

ydet|27rC| 


X 

Y 




( 8 ) 


where Co = (X • X'^), = (X • and 



All these quantities depend on the separation vector r only 
because of statistical homogeneity. This PDF is sufficient to 
compute the expectation of any quantity involving deriva¬ 
tives of the potential and the density up to second order. 
All the coefficients can easily be computed from the power 
spectrum of the potential 


27r 

d^ dk Pk (k) exp(z/c r cos 0 ) 
Uil+i2+jl+j2 + 'i- 

/i+i2(_j)ji+j2(cos6))'i+J'i(sin6»)'"+^'" - -, (9) 

and 

27r 

de dfcPfe(fc)fc^”+\ 

where the power spectrum of the potential Pk{k) can include 
a filter function on a given scale. In this work, we use a 
Gaussian filter defined in Fourier space by 

■ (i«) 

For instance, the one-point covariance matrix for 


{A4>,4>+,4>^,Ax, 

X^,X^) 

at a given point simply reads 


/ 1 

0 

0 

-7 

0 

0 \ 


0 

1/8 

0 

0 

- 7/8 

0 

C02 = 

0 

-7 

0 

0 

1/8 

0 

0 

1 

0 

0 

1 

0 

00 


0 

- 7/8 

0 

0 

1/8 

0 


V 0 

0 

- 7/8 

0 

0 

1/8 / 


where 7 = y^(n + 2)/(n + 4) for a scale-invariant density 
power-spectrum with spectral index n (i.e n — 4 for the po¬ 
tential). Note that the first derivatives of the density and the 
potential fields are decorrelated from the second derivatives 
meaning that 





Co = 


Coi 

0 


0 

C02 



Figure 5. Two-point correlation functions as a function of the 
separation r in units of the smoothing length for a power-law 
2 D power spectrum with spectral index n = — 1/2 i.e 7 = -\/ 3 / 7 . 
For aesthetic purpose, these functions have been rescaled by their 
value in r = 0 as explicitly written in the legend. 


of the potential and the density fields (0i, 02 , iC 2 ) as a 
function of 7 ' = cr|/cri /(73 



1/2 

0 

-772 

0 

\ 

0 

1/2 

0 

-772 



-772 

0 

1/2 

0 



0 

-772 

0 

1/2 

/ 


The two-point covariance matrix, can be similarly 
derived. In particular, its restriction to the second deriva¬ 
tives of the density and the potential fields can be written 
as a function of the nine ^ functions defined in equation (5) 
(for fg — (pcj), (j)x^ xx) (see Figure 5 and Appendix G 2 ): 



tA+ 

aAA 

^77 _ tX X 

4 

0 

tAA 

S^cc 

^A+ 

^(bx 

tAA 

_ tXX 

4 ^ 4 )X 

0 

\ 

fA+ 

0 

tA+ 

0 


0 

0 

X 

0 

0 

^X X 
^ 4 )x 


^AA 

tA+ 

^d)X 

0 

^AA 

^xx 

£A + 

Sxx 

0 


^A+ 

aAA 

^cj)x tXX 

4 ^(j)x 

0 

£A+ 

Sxcc 

tAA 

^xx tXX 

4 Scccc 

0 


0 

0 

^X X 
^ 4 >x 

0 

0 

^X X 
Scccc 

/ 


Once the joint PDF given by Equation ( 8 ) is known, 
it is straightforward to compute conditional PDFs (in par¬ 
ticular subject to a critical point constraint C(crit) = 

\ det{xij) \S-D{xi)). Given the conditionals, simple algebra 
then yield the conditional density and spin. More specifi¬ 
cally, relying on Bayes theorem, the conditional can be ex¬ 
pressed in terms of the joint PDF - equation ( 8 )- as 


P(X|Y,pk) 


7^(X,Y.pk) 

P(Y,pk) 


where 


V{Y,pk) = J dYP(X,Y)C(pk) 


where Coi is the one-pt covariance matrix of the gradients 


is the marginal distribution describing the likelihood of 
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a given peak, pk (the transverse cross-section of an infi¬ 
nite filament) with a given geometry. Once the conditional, 
7^(X|Y,pk) is known, it is straightforward^ to compute the 
expectation of any function, /(X) as 

(/(x)ipk>=ydxp(xipk)/(x), (11) 

which, when /(X) is multinomial in the components of X 
can be carried out analytically. In the following, we will con¬ 
sider in turn functions which are indeed algebraic function 
of X. 


3.1.2 Constrained density maps 


From Equation (8), given a contrast u and a geometry for 
the saddle (or any critical point) defined by = Ai — A 2 , /i = 
A 1 -I-A 2 (where Ai > A 2 are the two eigenvalues of the Hessian 
of the density field H - both negative for a peak), the mean 
density contrast, (5ext = ((5|ext), (in units of ( 72 ) around the 
corresponding critical point can be analytically computed 


(5(r|ext, K, /i, jy) 


1 — 7 ^ 

+ 4(fT-H.f)ct+, (12) 


where H is the detraced Hessian of the density and f = r/r 
so that 


f H f = 


cos(2^) 

2 


(13) 


r being the distance to the critical point and 0 the angle from 
the eigen-direction corresponding to the first eigenvalue Ai 
of the critical point. When r goes to zero, given the prop¬ 
erties of the ^ functions (see Figure 5), the density trivially 
converges to the constraint ly. 


3.1.3 Constrained 2D spin maps 
In two dimensions, the spin is a scalar given by 
'S2;(j(’) — ^ ^ 

where e = Cijs is built upon the totally anti-symmetric rank 
3 Levi-Civita tensor Cijk- Since equation (14) is quadratic in 
the fields x and 0, equation (11) can be readily applied to 
compute analytically its conditional expectation. The angu¬ 
lar momentum generated by TTT as a function of the polar 
position (r, 0) subject to the same critical point constraint at 
the origin with contrast z/, and principal curvatures (Ai, A 2 ) 
is given by the sum of a quadrupole ((x sin 2^) and an oc- 
tupole (oc sin 4^) 

(s^|ext) = s^(r|ext,z^,/i,z/), (15) 

= -16(f'^- e • H • f) (7^ + 2(f'^- H • f)s7), 

where the octupolar coefficient can be written as 



^ see http://tinyurl.com/inmbse3z which describes an imple¬ 
mentation in MATHEMATICA of the conditional probability. 


and the quadrupolar coefficient s^P reads 


+ 


while 


u - sin(26') 

r • e • H • r = —K-^^ 


( 16 ) 


Equation (15) is remarkably simple. As expected, the spin, 
s^, is identically null if the filament is axially symmetric 
(/Y = 0). It is zero along the principal axis of the Hessian 
(where ^ = 0 mod 7r/2 for which f^-e-f = 0). Near the peak, 
the anti-symmetric, sin(2^), component dominates, and the 
spin distribution is quadrupolar. For scale-invariant density 
power spectra with index n (n — 4 for the potential), Sz can 
be computed explicitly. At small separation, Sz behaves like 


Sz OC /^((n + 2)iy + \/(n + 2)(n + 4)/i)r^ sin(2^), (17) 

which shows explicitly that the quadrupolar term domi¬ 
nates. Figure 6 displays the mean density and spin map for 
a power-law power spectrum with index n = 1/2 around a 
2D peak of the density field with geometry z/ = 1, Ai = — 1 
and A 2 = —2. 

At this stage it is interesting to understand how much 
angular momentum is contained into spheres of increasing 
radius that would feed the forming object at different stages 
of its evolution. For instance let us assume there is a small- 
scale overdensity at (one of the four) location of maximum 
angular momentum (denoted r* hereafter) and let us filter 
the spin field with a top-hat window function centered on r* 
and of radius Rth . The resulting amount of angular momen¬ 
tum as a function of this top-hat scale is displayed in Fig¬ 
ure 7. During the first stage of evolution, the central object 
will acquire spin constructively until it reaches a Lagrangian 
size of radius T^th = r** and feels the 2 neighbouring quad¬ 
rants of opposite spin direction. The spin amplitude then 
decreases and becomes even negative before it is fed by the 
last quadrant of positive spin. The minimum is reached for 
radius around 2.4r*. This result does not change much with 
the contrast and the geometry of the peak constraint. Fig¬ 
ure 7 is the Lagrangian counterpart of Figure 4 of Laigle 
et al. (2015) (or Figure 7 of Pichon & Bernardeau (1999)) 
which displays the quadrant of vorticity in the vicinity of 
filaments. 


3 . 1.4 Cosmie varianee on spin 

On top of the mean spin, one can also compute the disper¬ 
sion of the spin described by 

a(r) = /(si(r)) - (s,(r))^ (18) 

A map of this spin dispersion is shown on Figure 8. Compar¬ 
ing Figure 8 to Figure 6, we see that spin direction fluctuates 
along the major axis of the filament cross section, and best 
defined along its minor axis. As the conditional statistics is 
Gaussian, the whole spin statistics (third moments, ...) can 
in principle be similarly computed. 
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Figure 6. Left: mean density (contrast) field near a 2D peak of height u = 1, Ai = —1 and A 2 = —2 for a power spectrum with index 
n = II 2 computed from Equation (12). Contours are displayed from (5 = —1 to 1 by step of 1/4 as labeled. The x and y axes are 
in units of the smoothing length. Right: corresponding mean spin colour coded from blue (negative) to red (positive) computed from 
Equation (14). The flattening of the filament’s cross section induces a clear quadrupolar spin distribution in its vicinity. 


L 



Figure 7. Evolution of the amount of algebraic angular momen¬ 
tum in sphere of radius Rth centered on r* The density power 
spectrum index is n = —3/2, the height of the peak in (0,0) is 
u = 1 and principal curvatures Ai = — 1, A 2 = —2. The amplitude 
of the spin is normalised by its maximum value around Rth = 


3.1.5 Zel’dovich mapping of the Spin 

Figure 9 displays the image of the initial density field (resp. 
initial spin field) translated by a Zel’dovich displacement. 
The displacement is proportional to (0i,02) here and its 
expectation given a central peak is trivially computed from 
the conditional PDFs. The resulting quadrupolar caustics is 
qualitatively similar to the quadrupolar geometry of the vor- 
ticity field measured in numerical simulations (Laigle et al. 
2015). Indeed, as discussed in that paper, there is a dual 
relationship between such Eulerian vorticity maps and the 
geometry of the spin distribution within the neighbouring 
patch of a 3D saddle point. 



Figure 8. 2D spin dispersion (defined in Equation (18)) near a 
2D peak of height n = 1 and curvatures Ai = — 1 and A 2 = —2 
for a power spectrum with index n = 1/2 


3.2 Transition mass for long filaments 

Up to know we assumed that the geometry of the critical 
point was given. Let us now build the joint statistics of the 
spin and the mass near 2D peaks. 
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X / Rs 
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-2 


- 2-10 1 2 
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Figure 9. Left: Stream lines of the 2D velocity field (defined as the potential gradient) near a 2D peak of height 1 ^ = 1 and curvatures 
Ai = — 1 and A 2 = —2 for a power spectrum with index n = ll2. Right: Zel’dovich mapping of the spin distribution. There is a good 
qualitative agreement between the vorticity section presented in Laigle et al. (2015) and this spin map. 


S.2.1 Geometry of the most likely eross seetion 


same orientation, can then simply be expressed as 


Let us now study what should be the typical geometry of a 
peak. Following Pogosyan et al. (2009), it is straightforward 
to derive the PDF for a point to have height u and geometry 
as in their notation J2 = so that 


V{v,K,h) 


yi -72 



V + 7/1 

71 -P 




Now the PDF for a yeak to have height v and geometry 1\ 
becomes: 


^ = Trr* , (20) 

where r* = r*(i/, k) is the position of a maximum of angular 
momentum from the peak. Because of the quadrupolar anti¬ 
symmetric geometry of the angular momentum distribution 
near the saddle point, it is typically twice as small (in units 
of the smoothing length) as one would naively expect. 

For power-law density power spectrum with spectral in¬ 
dex in the range n g] — 2 , 2], a good fit to its scaling is given 
by 


V{iy, Hi,Ii\pk) 


V3k\(Ii - k){Ii + k)\ ^ 

27r^l — 7 ^ 


exp 


1 f y + 7/1 

2 (71 - 7^ 




The maximum of this PDF is trivially reached for 
R = y/l/3 and h — —y/7/3. 


(19) 
R — 


3.2.2 The size and area of eonstant polarity quadrants 

From equation (15), it appears clearly that the extension of 
the region of influence of the critical point is limited, and 
peaks within each quadrant at some specific (r*, ^*) position. 
Moreover, for small enough Av, the quadrupole dominates, 
and the extremum is along ^ = 7r/4. It is therefore possible 
to use r* to define an area in which the spin is significantly 
non zero within each quadrant. Let us compute r*, as the 
radius for which Sz{0 = 7r/4) is maximal as a function of . 
The area of a typical quadrant, in which the spin has the 


^ Setting 6 = 7T14 effectively neglect the octupolar part of Sz- 


Rs 


3 

250 


(n-5)^ 


+ 


10 ’ 


( 21 ) 


where r* was computed for the mean geometry given by 
V = ./fji'y, K = 7V3 and h = -s/Pp>- 


3.2.3 Critieal mass sealing 

The critical mass is the mass of maximum spin alignment. 
In simulation, it has been shown by Laigle et al. (2015) to 
be Merit ~ IO^^Mq at redshift 0. The authors claimed that 
the critical mass is related to the mass contained in a typi¬ 
cal quadrant of vorticity. In this work, we have computed in 
Lagrangian space the typical area of a quadrant (see Equa¬ 
tion 20). This area is a function of the smoothing scale. In 
order to compute it, we need to define a scale. It is reason¬ 
able that the maximum spin alignment should be reached 
for filament that has just collapsed at redshift 0. Indeed, for 
larger scale filaments, part of the haloes do not lie inside the 
filament but in the nearby wall which will therefore decrease 
the mean spin-filament alignment. In previous sections, we 
focused on u = 0.9 filaments. The model of the cylindrical 
collapse then say that those filaments have just collapsed at 
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redshift 0 for a top-hat initial smoothing scale cr(i? th) = 1.6 
which corresponds to a smoothing length Rtu — 2.2Mpc//i. 
We can therefore compute the corresponding r* which is 
r* ~ l.QRs ~ O.Ti^TH = 1.5Mpc//i and corresponds to the 
mass Merit = ~ 1.510 ^^Mq, in good agreement 

with the value measured in simulations. Its redshift evolu¬ 
tion is also predicted by the formalism through the cylindri¬ 
cal collapse and could be compared to simulation in future 
works. 

Note that this line of reasoning could be made more 
rigorous by adding new ingredients in the formalism: a peak 
constraint at the location of the spin with a smoothing 
length Rh < Rs so that one can vary Rh (without any as¬ 
sumption on the additivity property of the spin) and see 
how the spin changes. This formalism can be implemented 
in two dimensions (see Appendix A2 and A3) and leads to 
the same order of magnitude for r*. 


4 3D SPIN NEAR AND ALONG FILAMENTS 

Let us now turn to the truly three dimensional theory of tidal 
torques in the vicinity of a typical filament saddle point. Be¬ 
yond the obvious increased realism, the main motivation is 
that the 3D saddle theory fully captures the mass transition. 

In three dimensions, we must consider two competing 
processes. If we vary the radius corresponding to the La- 
grangian patch centered on the running point, we have a 
spin-up (along e^) arising from the running to wall-running 
to saddle tidal misalignment and a second spin-up (along 
G(f)) arising from running to filament- running peak tidal 
misalignment. To each position in the vicinity of the central 
saddle point, we can assign M(r) together with cos /iz(r) and 
cos/i0(r), the cosines of the angle between the spin of the 
patch and the and direction respectively. Eliminating 
r yields cos/A z{M) and cos and therefore yields an 
estimate of the transition mass. 


4.1 Spin distribution along and near filaments 

The formalism developed in Section 3 can easily be extended 
to three dimensions. A critical (saddle) point constraint is 
now imposed. This critical point is defined by its geometry, 
namely its height v and eigenvalues Ai ^ A2 ^ A3. Note 
that such a critical point is a filament-type saddle point if 
Ai ^ 0 ^ A2 ^ A3. In what follows, we decouple the trace 
from the detraced part of the density Hessian and therefore 
define the three curvature parameters /i = Ai + A2 + A3, 
= Ai — A2 and /^2 = A2 — A3. 


4 . 1.1 Mean density field around a eritieal point 

The resulting mean density (contrast) field subject to that 
critical point constraint becomes (in units ol a 2 ) : 


(5(r|crit, /i, a^i, /^ 2, h>) 


1 — 7 ^ 1 — 7 ^ 

+ y (f^-H-f)7+, (22) 



Figure 10. Two-point correlation functions as a function of the 
separation r in units of the smoothing length for a power-law 3D 
power spectrum with spectral index n = — 2 i.e 7 = \/3/3. As 
in two dimensions, the correlations are rescaled at the origin. See 
also Figure Cl 


where again H is the detraced Hessian of the density and 
f = r/r and we define in 3D as 

7/ = (Ax0+), (23) 

with 0+ = 011 — (022 + 033)/2. The other ^ functions are 
defined in the same way as in two dimensions (see equa¬ 
tions 5) and displayed on Figure 10. Note also that r"*" • H • r 
is a scalar defined explicitly as ViHijTj. Figure 11 dis¬ 

plays the mean density field around a typical filament-type 
saddle point. The elongation of the filament along the Oz 
axis together with the flattening of the filament in the plane 
of the wall (Oxz) are clearly-visible on this figure. 


4 . 1.2 Mean spin field around a eritieal point 

As in two dimensions, the expected spin can also be com¬ 
puted. In three dimensions, the spin, s, is a vector which 
components are given by 

Si — ^ ^ CijkXklf^lj 5 (^^) 

j,k,l 

with e the rank 3 Levi-Civita tensor. It is found to be or¬ 
thogonal to the separation and can be written as the sum of 
two terms 


s(r|crit, ia) = —15(s^^^ + s^^^) • (r"*"- e • H • r), (25) 

where is a scalar operator that depends on the height u 
and trace of the Hessian R 



+Y^ , 
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Figure 11. Mean density around a filament saddle point of height u = 1.257, = 0.31, A 2 = —0.56, A 3 = —1 for a power-law 3D 

power spectrum with spectral index n = —2 i.e 7 = \/3/3. The Ox^ Oy^ and Oz axes are in units of the smoothing length and z is the 
direction of the filament while the wall is in the plane Oxz. Contours represent the isocontours of the density S = —0.3, 0 , 0.5 from blue 
to red. A 3D view is displayed on the left panel and a cut in the plane Oyz is shown on the right panel. 


and a combination of a matricial and a scalar operator 
that depends on the detraced part of the Hessian 

=-I [2(gt^ - 

+ + 5^7+)47) (f^ • H • f)l3] 

with I3 the identity matrix, operating on the vector 

r^- e • H • r = ^ ridjk'H.kih ■ (26) 

ikl 

Note that the dependence with the distance r is encoded in 
the two-point correlation functions, while the geometry 
of the critical point is encoded in the terms corresponding 
to the peak height, trace and detraced part of the Hessian 
and the orientation of the separation is in f. Equation (25) 
is also remarkably simple: as expected the symmetry of the 
model induces zero spin along the principal directions of 
the Hessian (where r"*"- e • H • r = 0) and a point reflection 
symmetry (f —f). Note that the correlation functions, ^ 
can be evaluated for arbitrary power spectra (such as power 
laws, see Appendix C2, or ACDM, see Appendix C3), hence 
Equation (25) is completely general. 

Eor scale-invariant density power spectra with index n 
(n — 4 for the potential), s can be computed explicitly. At 
small separation, the term proportional to • H • r goes like 
and is thus negligible compare to the rest (that scales 
like r^). The spin coordinates in the frame of the Hessian 
are therefore quadrupolar 

s oc IP, n)yz, g{Xi, n, n)xz, h{Xi, u, n)xy) . (27) 

Eigure 12 illustrate the mean spin geometry around 
a typical saddle point. All the symmetry properties (anti¬ 
symmetry, octopole, ...) described in this section are clearly 


seen on this figure. In the plane of the saddle point, spins are 
aligned with the filament direction. When moving towards 
the nodes, the spins become more and more perpendicular 
(and more and more along e(^). 


4 . 1.3 Cosmic variance on 3D spin 

It is of interest to also study the variance of the spin align¬ 
ment cr(r|ext, /i, /^i, /^2, k') defined as 

a — \J (cos^ 0 ) — (cos^)^ , (28) 

where cos 0 = s • Gz /11 s 11. It requires the numerical evalua¬ 
tion of a 12D integral. In contrast, the mean of the spin s (as 
computed in Section 4.1.2) or its square can be analyt¬ 
ically computed. We therefore propose to approximate the 
dispersion of the spin alignment with the following related 
estimator 


a = 


(^i) - (^4^ 

(s-s> 


(29) 


where Sz is the component of the spin along the z-axis i.e 
along the filament direction. Eor the sake of readability, we 
do not write down the result of the integration here but 
display in Eigure 13 the map of the alignment dispersion 
a around a typical saddle point. This standard deviation is 
roughly constant around 0.6 and decreases to 0.3 in the 
close vicinity of the saddle point. Note that the spin direction 
is again best defined along its minor axis. This would be the 
best place to measure spin alignments in observations. 
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Figure 12. The velocity and Spin flow near a vertical filament (in 
red) embedded in a (purple) wall for (x, y, z) G [—2i?s, 2Rs]^. The 
purple and green flow lines trace the (Lagrangian) 3D velocities 
(upwards and downwards respectively). The red and blue arrows 
show the spin 3D distribution, while the three horizontal cross 
sections show spin flow lines in the corresponding plane. Note 
that the spin is along Sz in the mid plane and along away 
from it, and that it rotates in opposite direction above and below 
the mid-plane. See also the interactive version at http://www. 
iap.fr/users/pichon/AM-near-saddle.html 


4.2 Mean saddle point geometry 

Here we want to compute the mean values of i/, Ai < A 2 < 
0 < A 3 of a typical saddle-point of filament type. Let us 
start from the so-called Doroshkevich formula for the PDF 
of these variables: 


(Ai) ~ —0.87, (A 2 ) ~ —0.40 and (A 3 ) ~ 0.51. However, this 
saddle point does not belong to the skeleton of the density 
field but to its inter-skeleton (see Pogosyan et al. (2009)). 
We thus want to impose an additional constraint which 
is A 2 + A 3 < 0. Let us call those saddle points “skeleton 
saddles”. The PDF at those points becomes 


V{iA, Ai|skl) 


26460y^P(!^,Ai)/3e(A3)^, , , , 

1421x/2 - 735V3 + 66n/42®^ ' 


(31) 


The expected value of the density and the eigenvalues at a 
skeleton saddle position now becomes (ly) ~ 1.257, (^ 1 ) ~ 
- 1 . 0 , (A2) ^ -0.56 and (A3) 0.31. 


4.3 Spin flip : from spatial to mass transition 

The geometry of the spin distribution near a typical skele¬ 
ton saddle point (as defined by equation (31)) allows us to 
compute the mean alignment angle between the spin and 
the filament (see Section 4.3.1 below). In turn, the shape of 
the density profile in the vicinity of the same critical point, 
together with an extension of the Press-Schechter theory in¬ 
volving a filament background split, allows us to estimate 
the “typical” mass of the dark matter haloes forming in any 
spatial position around the saddle point (Section 4.3.2 be¬ 
low). The alignment-angle map and the typical-mass map 
will together yield a prediction for the transition mass. 


4 . 3.1 Spin flip along filaments 

Section 4.1 showed that the mean spin flips from alignment 
in the plane of the saddle point to orthogonality when going 
towards the nodes. This can be quantihed by measuring the 
curvilinear coordinate along the hlament at which the spin 
flips. 

Let us consider the mean modulus of the projection of 
the spin along the Gz and axes within a plane of height 
z 

^z/<i,{A = J d®dj/|s^/^(r)|/||s(r)||. (32) 

Figure 14 displays Szjcpiz) as a function of along the hla¬ 
ment. Let us dehne 0 the flip angle so that 


,, 135(5/27r)"/' 


4v^l - 72 




X (A 3 — Ai)(A 3 — A 2 )(A 2 — Al) , 


where ( = (jy 7/i)/y^r^^^, /i = Ai -f A2 + A3, I2 = 
A1A2 + A2A3 + A1A3 and I3 = A1A2A3. Subject to a saddle- 
point constraint, this PDF becomes 

P(..A.|,ad) , ^”7^;;;^> /,9(A.)e(-A.). (30) 

after imposing the condition of saddle point 
I det ft^j(5|(5z:)(V(5)©(A3 )©(—A 2 ) for which as the gra¬ 
dient is decoupled from the density and the Hessian, 
only the condition on the sign of the eigenvalues and the 
determinant contribute. From this PDF, it is straightfor¬ 
ward to compute the expected value of the density and 
the eigenvalues at a saddle point position: {iy) ^ O.767, 


cos^(2:) = — ^z{z) _ 1 ^ 

In Figure 14, this flip angle is found to occurs around 
z = 1.5i?s which is very close to the r* measured in two 
dimensions (see Section 3.2). 

Alternatively, one can also compute at each position the 
mean alignement with the hlament direction Gz 


cos ^(r) 


s(r|crit) • Gz 
|| s(r|crit)|| 


(34) 


The result is shown on the right panel of Figure 15. Spins 
tend to align with the hlament (region in red) in the plane 
of the saddle point and becomes perpendicular to it when 
moving towards the nodes (region in blue). This is a transi¬ 
tion in Lagrangian space. Section 4.3.2 shows how to convert 
it into a transition in mass. 


MNRAS 000, 000-000 (0000) 
















14 


S. Codis, C. Pichon and D. Pogosyan 






0 

X / Rs 


0 

X / Rs 


Figure 13. Alignment dispersion a (defined in Equation 29) around a typical filament-type saddle point of height u = 1.257, = 

0.31, A 2 = —0.56, A 3 = —1 for a power-law 3D power spectrum with spectral index n = —2. The left panel displays a cut along the 
plane, Oxy, of the saddle point and right panel along the plane, Oxz, of the wall. The uncertainty on spin direction is smallest near the 
saddle. 



Figure 14. measure the height 2 ; corresponding to the transition 
from aligned to perpendicular to the filaments. The amplitude in 
both direction is averaged by plane following equation (32). The 
transition curvilinear coordinate is ztr = l.5Rs. 


4 . 3.2 Halo mass gradient along filament 

The local mass distribution of haloes is expected to vary 
along the large-scale filament due to changes in the underly¬ 
ing long-wave density. In the linear regime, the typical over¬ 
density near the end points (nodes) of the filament, where it 
joins the protocluster regions, may exceed the typical over¬ 
density near the saddle point by a factor of two (Pogosyan 
et al. 1998). During epochs before the whole filamentary 
structure has collapsed, this leads to a shift in the hierarchy 
of the forming haloes towards larger masses near the filament 
end points (the clusters) relative to the filament middle point 
(the saddle). This can be easily understood using the for¬ 
malism of barrier crossing (Peacock & Heavens 1990; Bond 


et al. 1991; Paranjape et al. 2012; Musso & Sheth 2012), 
which associates the density of objects of a given mass to 
the statistics for the random walk of halo density as the field 
is smoothed with decreasing filter sizes. Specifically these 
authors predict the first upcrossing probability for the crit¬ 
ical threshold at the filter scale corresponding to the mass 
of interest. The precise outcome of the formalism depends 
on the spectral properties of the field and the form of the 
smoothing filter, however it is clear that, in general, decreas¬ 
ing the barrier threshold increases the probability that such 
first upcrossing will happen at large smoothings, i.e large 
mass. A larger fraction of the Lagrangian space will then 
belong to large-mass haloes, at the expense of the low-mass 
ones. 

Following the presentation of Paranjape et al. (2012) 
of the Peacock-Heavens (Peacock & Heavens 1990) approx¬ 
imation -that was found to fit numerical simulations rather 
well, the number density of dark haloes in the interval 
[M, M + dM] is 

where f{a^,Sc) is given by the function 
/(c^^(5c) = exp 

f 2 dp(a^,Sc) 1 / 2 r Ni / 2 r A 

X i-cr - ^ - -p{a ,Sc)\np{a , ( 5 c) 1 . 

Here is the variance of the density fluctuations 
smoothed at the scale corresponding to M and p{a‘^,Sc) = 
1/2 (1 + erf((5c/\/2cr)) is the probability of a Gaussian pro¬ 
cess with variance to yield value below some critical 
threshold Sc. In Equation (36), P is the parameter depen- 


^lnp{s',Sc) ) X 


(36) 
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X / Rs X / Rs 

Figure 15. Left: cross section of Mp(r,z) (in units of IO^^Mq) along the most likely filament and in the direction x = y. Right: 
corresponding cross section of {cos 6 ){r,z), the normalised component of the spin aligned with the filament. The black dot represents 
the position of the saddle point. The mass of haloes increases towards the nodes, while the spin flips. 


dent on the filtering scale and, to less extend the underlying 
power spectrum, that specifies how correlated the density 
values at the same point when smoothed at different scales 
are. For Gaussian filter, the value F 4 is advocated. 

The overall mass distribution of haloes is well described 
by the choice dc = 3/5 (37r/2)^^^ = 1.681, motivated by 
the so-called spherical collapse model. When haloes form 
on top of a large-scale structure background, however, the 
long-wave over-density S{z) adds to the over-density in the 
proto-halo peaks. The effect on halo mass distribution, in 
this so-called peak-background split approach, can be ap¬ 
proximated as a shifted threshold Sc{z) = 1.681 — S{z) for 
halo formation. In Figure 16, we show that, as expected, 
the result of this long wavelength mode is a shift of the halo 
mass distribution towards larger masses. This shift can be 
characterized by the dependence on the threshold of M*((5c), 
defined as = Sc, or of the mass Mp{Sc) that corre¬ 

sponds to the peak of f{a‘^,Sc), i.e. the variance crp{z) de¬ 
fined by 

(Tp(z) = argmax(/((7^, (5o(2))^ • (37) 

Figure 17, right axis, shows these two characteristic vari¬ 
ances as functions of the threshold. Sc- 

The link to cosmology is established by relating the vari¬ 
ance cr^ to the mass of the objects. If a background held is 
absent, the variance is just the integral of the power spec¬ 
trum P{k) smoothed over a sphere of mass M 

cr^ = a^{M,Z) = D^{Z) J k‘^dkPik)WTH ((3M/47rp)^/®) 

where D(Z) is the linear growing mode of perturbations 
as a function of redshift Z and Wtu is the top-hat hlter. 
However, when large-scale structures are considered as hxed 
background, the variance of the relevant small-scale density 





Figure 16. Mass distribution for three values of density thresh¬ 
old, Sc = 1.681,1,0.681 from left (yellow dashed line) to right 
(blue solid line). The displayed function /(cr^) is defined in Equa¬ 
tion 36. 

huctuations that are responsible for object formation is re¬ 
duced, approximately as 

-a"(MLss), (38) 

where (T^(Mlss), given as well as a^{M) by Equation (38), 
is the unconstrained variance at the scales at which we 
have dehned the background large-scale density. This cor¬ 
rection is negligible when there is distinct scale separation 
between non-linear forming objects and the large-scale den¬ 
sity, i.e S{x) <C 1.681 but becomes important, truncating the 
mass hierarchy at Mlss, whenever large-scale structures are 
themselves non-linear. 

On Figure 17, left axis, the variances are converted into 
masses, M* and Mp according to equation (38). We choose 
here erg = 0.8, Z = 0, we dehne the mass in a Mpe 
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Figure 17. Characteristic variances (plain) and M* (yellow 
dashed) and Mp halo masses (blue dashed) as functions of the 
large-scale density in the peak-background split picture. 


comoving sphere for the best-fit cosmological mass density 
and we approximate the spectrum with a power-law of in¬ 
dex n = — 2, which allows to solve Equation (38) explicitly, 
giving the M(cr) relation as 

M{a,Z) = 2.6 X lO^^Me • (39) 

We consider filaments to be defined with R = 5/i“^Mpc 
Gaussian smoothing, which gives cr^(MLss) ~ 0.66. The evo¬ 
lution of Mp(S,Z) follows from putting equation (37) into 
equation (39). 


4-3.3 Spin orientation versus mass 

From the above described Mp-S relation, one can attribute 
a mass to each position depending on the value of the mean 
density at that location. The result is illustrated in Figure 15 
where the left and right panels display respectively the mass 
map and the spin alignment map around a typical saddle 
point. Eliminating the spatial position, r, between these two 
maps yields (cosO) as a function of Mp as shown on Fig¬ 
ure 18. The transition mass, Mtr for spin flip ((cos 0) = 0.5) 
is found to be of the order of 4 lO^^M©, assuming a smooth¬ 
ing scale of 5 Mpc//i, as used in Codis et al. (2012). This 
mass is in qualitative agreement with the transition mass 
found in that paper, all the more so as the redshift evolu¬ 
tion of this transition mass will also be consistent (scaling 
as the mass of non-linearity). 

It is quite striking that the geometry of the saddle point 
alone allows us to predict this mass. The two main ingre¬ 
dients for success are the point reflection symmetry of the 
spin distribution near the most likely filament-like saddle 
point on the one hand, and the peak background split mass 
distribution gradient along the filament towards the nodes 
of the cosmic web on the other hand. 


<cos 6 > 



Figure 18. Mean alignment as a function of mass for a smoothing 
scale for filaments of 5 Mpc/h. Error bars represent the error on 
the mean cosine in each bin of mass for the region [—2i?s, 2Rs] x 
[—2Rs, 2Rs] X [—2Rs, 2Rs] around a typical filament saddle point. 
The flip transition mass corresponds to 4 10^^ Mq. 

5 STATISTICS 

Up to now, we have considered the neighbourhood of a given 
unique typical saddle point as a proxy for the behaviour 
within a Gaussian random field (GRF). In view of our find¬ 
ing let us now first analyze the statistics of alignment for 
GRF, and then for fields corresponding to their simulated 
cosmic evolution down to redshift zero. 

5.1 Validation on GRF 

Let us consider the following experiment. Let us generate 2D 
or 3D realizations of GRF smoothed on two successive scales, 
Lh and Ls ^ Lh- In the first maps, let us build a catalogue of 
positions, Vh and heights, Uh corresponding to “small-scale” 
peaks. From the second maps, let us identify the loci, of 
the corresponding “large-scale” peaks (in 2D) and (filament 
type) saddles (in 3D), and build the corresponding fields s(r) 
(via fft using equation (24)). This field allows us to assign 
a spin to each ‘halo’ at position Vh and a closest saddle, r^. 
Given the relative position Vh — Ts as measured in the frame 
defined by the Hessian at r^, we may project the direction 
of the spins, s = s/s of all ‘haloes’ in the vicinity along the 
corresponding local cylindrical coordinate (ei^, e^). We 

may then compute the one point statistics of = s • per 
octant. 

5.1.1 2D GRF fields spin flip 

In two dimensions, the expectation is that the spin should be 
aligned or anti-aligned with Cz depending on each quadrant. 

Let us first start with a set of 25 2D 2048^ maps from a 
power-spectrum with n = —1/2. The map is first smoothed 
with Gaussian filter of width Lh — ^ pixels, and the posi¬ 
tions of the peaks are identified. It is then smoothed again 
over Ls — 24 pixels, exponentiated (in order to mimick 
the almost log-normal statistics of the evolved cosmic den¬ 
sity field), and the corresponding Hessian and tidal fields 
are computed, together with the momentum map, which is 
thresholded above 1/30^^ of its highest value (see Figure 19). 
The peaks of this second map are identified as ’saddles’ for 
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Figure 19. Top: example of spin map (colour coded by sign) 
generated following the prescription of Section 5.1.1. Bottom: the 
local frame (in red, long axis and blue) around a couple of ‘sad¬ 
dle’s. The black contours correspond to the density. 


contrasts higher than 2.5. Figure 20 shows that the average 
spin of ‘haloes’ in each quadrant is flipping from one quad¬ 
rant to the next, with a statistically significant non zero 
mean value in each quadrant. 

5.1.2 3D GRF fields spin flip 

Let us similarly consider a set of 20 three dimensionnal 256^ 
cubes from a power-spectrum with n = — 2. The cube is first 
smoothed with Gaussian filter of width = 4 pixels, and 
the positions of the peaks are identified. It is then smoothed 
again over Lg = 24 pixels, exponentiated, and the corre¬ 
sponding Hessian and tidal fields are computed, together 
with the spin field, which is thresholded above 1/30^^ of its 
highest value. The saddle of this second cube are identified 
as for contrasts higher than 1. Only peaks closer than one 
smoothing length from the large-scale saddles are kept. The 
angle between their spin and the filament axis is computed 
and stored depending on the octant they belong to. In this 
section, the octants are numbered from 1 to 8 depending 
on the separation from the peak to the saddle r = {x,y,z): 
x,y,z > 0 (#1), x<0&7/, z>0 (#2), x,y < 0 L z > 0 
(#3), y < 0 k x,z > 0 (#4), z < 0 k x,y > 0 (#5), 
x,z<0 &:y>0 (#6), x,y,z < 0 (#7) and y,z<0Szx>0 



Figure 20. Alignment of ‘spin’ along in two dimensions as 
a function of quadrant rank, clockwise. As expected, from one 
quadrant to the next, the spin is on average unambiguously flip¬ 
ping sign. 



Figure 21. Alignment of the spin along the filamentary direction 
depending on the considered octant. As predicted by the theory, 
the z-component of the spin is flipping sign from one octant to 
the other. 


(#8). Fig. 21 shows that, as expected, the component of the 
spin aligned with the filament axis is flipping sign from one 
octant to the other. 


5.2 Validation on dark matter simulations at z = 0 

Let us now identify the Eulerian implication at redshift zero 
of the above sketched Lagrangian theory. For this we must 
rely on N-body simulations. Hence we now make use of the 
43 million dark matter haloes detected at redshift zero in 
the Horizon 47r N-body simulation (Teyssier et al. 2009) to 
test some of the outcomes of the Anisotropic Tidal Torque 
Theory presented in this paper. This simulation contains 
4096^ DM particles distributed in a 2 h~^Gpc periodic box 
and is characterized by the following ACDM cosmology: 
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fim = 0.24, Qa = 0.76, n = 0.958, Ho = 73 km-s“^-Mpc“^ 
and as = 0.77 within one standard deviation of WMAP3 
results (Spergel et al. 2003). The initial conditions were 
evolved non-linearly down to redshift zero using the adap¬ 
tive mesh refinement code RAMSES (Teyssier 2002), on a 
4096^ grid. The motion of the particles was followed with 
a multi grid Particle-Mesh Poisson solver using a Cloud-In- 
Cell interpolation algorithm to assign these particles to the 
grid (the refinement strategy of 40 particles as a threshold 
for refinement allowed us to reach a constant physical reso¬ 
lution of 10 kpc, see the above mentioned two references). 

The Friend-of-Friend Algorithm (Huchra & Geller 1982) 
was used over 18^ overlapping subsets of the simulation with 
a linking length of 0.2 times the mean interparticular dis¬ 
tance to define dark matter haloes. In the present work, we 
only consider haloes with more than 40 particles (the parti¬ 
cle mass being 7.7 x 10^Mq). The mass dynamical range of 
this simulation spans about 5 decades. 

The filament’s direction is then defined via the global 
skeleton algorithm introduced by Sousbie et al. (2009) and 
based on Morse theory. It defines the skeleton as the set of 
critical lines joining the maxima of the density field through 
saddle points following the gradient. In practice Sousbie 
et al. (2009) define the peak and void patches of the den¬ 
sity field as the set of points converging to a specific local 
maximum/minimum while following the field lines in the di¬ 
rection/opposite direction of the gradient. The skeleton is 
then the set of intersection of the void patches i.e. the sub¬ 
set of critical lines connecting the saddle points and the local 
maxima of a density field and parallel to the gradient of the 
field. In practice, the ~70 billion particles of the Horizon-47r 
were sampled on a 2048^ cartesian grid and the density field 
was smoothed using mpsmooth (Prunet et al. 2008) over a 
scale of 5 h~^Mpc corresponding to a mass of 1.9 x 
This cube was then divided into 6^ overlapping sub-cubes 
and the skeleton was computed for each of these sub-cubes. 
It was then reconnected across the entire simulation volume 
to produce a catalog of segments which locally defines the 
direction of the filaments. 

Figure 22 demonstrates that the spins of the 43 million 
dark haloes of the simulation obey the expected mass de¬ 
pendent flip predicted by the theory presented in Section 4. 
On top of the alignment with the filament direction found 
e.g in Codis et al. (2012), haloes are shown to have a spin 
increasingly perpendicular to at low-mass (red) and up 
to the critical mass (~ 10^^M©), while high-mass haloes 
have a spin parallel to the direction. The transition from 
alignment to orthogonality occurs around Mtr — 5- 

Figure 23 shows that the spins tend to be more aligned 
with the filament axis when getting closer to the saddle 
point. The alignment decreases from cosO = 0.511 at r cp 20 
Mpc/h to cos^ = 0.506 at r < 1 Mpc/h. This qualitative 
trend is in full agreement with the anisotropic tidal torque 
theory picture presented in Section 4 for which on average, 
spins are aligned with the filament axis in the plane of the 
saddle point and become misaligned when going away from 
this saddle point. 

Figure 24 displays the occupancy of haloes along the 
filaments. It appears that the higher the mass, the more 
concentrated they are far from the saddles. This is in good 
agreement with the halo mass gradient along the filaments 
described in Section 4.3.2. 



COS 0 


Figure 22. Alignment of ‘spin’ along in the Horizon-47r sim¬ 
ulation. The normalised histogram of the cosine of the angle be¬ 
tween the spins and the closest filament’s direction is displayed. 
Deviations from the ^ = 0 uniform distribution are detected and 
depends on the dark matter halo mass. Haloes have a spin aligned 
with the Cff) direction on average at low-mass (red) and perpen¬ 
dicular to it at larger mass (blue). 



Figure 23. Alignment of the spins and the filaments in the 
Horizon-47r simulation as a function of the distance to the clos¬ 
est saddle point (from red - OMpc/h - to blue -10 Mpc/h -). 
The alignment decreases with the distance to the saddle point as 
predicted by the anisotropic tidal torque theory model. 


Overall, the above GRF experiments as well as the re¬ 
analysis of the Horizon 47r N-body simulation seem consis¬ 
tent with the prediction of the theory presented in Section 3 
and 4. While the former demonstrates that interferences 
from neighbouring saddles do not wash out the tide correla¬ 
tions, the latter suggests that on the scales probed by this 
experiment, this Eulerian measure still captures features of 
the underlying Lagrangian theory. 
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distance to saddle (Mpc/h) 


Figure 24. (Log-) Fraction of haloes of different mass (from red 
to blue in solar mass units) as a function of the distance to the 
saddle point in the Horizon-47r simulation. Low-mass haloes (red) 
lie almost uniformly along the filaments (with a small concentra¬ 
tion -not clearly seen in logarithmic units- around the saddles 
due to consumption when going towards the nodes) while high- 
mass haloes (blue) are more concentrated far from the saddles. 


6 CONCLUSIONS AND PERSPECTIVES 

TTT was revisited while focussing on an anisotropic peak 
background split in the vicinity of a saddle point. Such crit¬ 
ical point captures as a point process the geometry of a 
typical filament embedded in a given wall (Pogosyan et al. 
1998). The induced mis-alignment between the tidal tensor 
and the hessian of the density simply explains the surround¬ 
ing transverse and longitudinal point reflection-symmetric 
geometry of the spin distribution near filaments. 

This geometry of the spin field predicts in particular 
that less massive galaxies have their spin parallel to the 
filament, while more massive ones have their spin in the 
azimuthal direction. The corresponding transition mass fol¬ 
lows from this geometry together with its scaling with the 
mass of non linearity, in good agreement with measurements 
in simulations. 

The main findings of this paper are: i) galaxies form 
near filaments embedded in walls, and flow towards the 
nodes: this anisotropic environment produces the long wave 
modes on top of which galactic haloes pass the turnaround 
threshold; ii) a typical filament is elongated and flattened: 
as a point process, it is therefore best characterized by its 
triaxial saddle points; iii) the spin geometry is octupolar 
in the vicinity of the saddle point, displaying a point re¬ 
flection symmetry; iv) the mean spin field is parallel to the 
filament axis in the plane of the saddle point and becomes 
azimuthal away from it; v) the constrained tidal torque the¬ 
ory presented in this paper allows to accurately predict the 
transition mass of the spin-filament alignment measured in 
simulations, vi) this theory seems consistent with both GRF 
experiments and results from N-body simulations, vii) a dual 
theory describes spin alignments in voids (see Appendix B). 


6.1 Discussion 

One of the striking features of this anisotropic extension of 
TTT is the induced quadrupolar point-symmetric flattened 
geometry of the spin distribution near a saddle point, which 
effectively scales down by one order of magnitude the transi¬ 
tion mass away from the mass of non-linearity, in agreement 
with the measured scaling. The qualitative analysis derived 
from first principles in the vicinity of a given saddle point 
seems to hold when considering realizations of GRF, once 
proper account of the induced geometry near such points 
is taken care of. In effect, we have shown that the geome¬ 
try of the saddle point provides a natural ‘metric’ (the local 
frame as defined by the hessian at that saddle point) rela¬ 
tive to which we can study the dynamical evolution of dark 
haloes along filaments. It should allow us to study how galac¬ 
tic feeding (via helicoidal cold flow, see Dubois et al. 2014) 
should vary with curvilinear coordinate along the filament. 
It was indeed found in that paper using hydrodynamical 
simulations that such flows were reaching galaxies in the so- 
called circum-galactic medium with velocities roughly paral¬ 
lel the polar axis. Taken at face value, such findings suggest 
that the flow feeding galaxies has significant helicity during 
that phase. 

Another striking feature of this Lagrangian framework 
is that it captures naturally the arguably non-linear Eulerian 
process of spin flip via mergers. Recently, Laigle et al. (2015) 
showed that angular momentum generation of haloes could 
be captured in Eulerian space via the secondary advection 
of vorticity which the formation of the filament generates, 
whereas we show in this paper that it may also be described 
in Lagrangian space via the analysis of the anisotropic tides 
generated by the filament to be. No description is more fun¬ 
damental than the other but are the two (Eulerian versus 
Lagrangian) sides of the same coin. The mapping between 
the two descriptions requires a reversible time integrator, 
such as the Zel’dovich approximation, which clearly limits 
its temporal validity to weakly non-linear scales. Our proxy 
for the spin, equation (3), is an approximation which seems 
to quantitatively capture the relevant physics. It is remark¬ 
able that such an (admittedly approximate) straightforward 
extension of TTT captures what seems to be the driving pro¬ 
cess of spin orientation acquisition and its initial evolution. 
It is also sticking that very simple closed form for the spin 
orientation distribution in the vicinity of the saddle point 
are available for this proxy. 

Our theory here has focussed on a two-scale process. 
Given the characteristics of AGDM hierarchical clustering 
one can anticipate that this process occurs on several nested 
scales at various epochs - and arguably on various scales at 
the same epoch. The scenario we propose for the origin of 
this signal is, like the signal itself, relative to the linear scale 
involved in defining the filaments and as such, multi-scale. 
It will hold as long as filaments are well defined in order 
to drive the local cosmic flow. In other words, one expects 
smaller-scale filaments are themselves embedded in larger- 
scale walls. The induced multi-scale anisotropic flow tran¬ 
spires in the scaling of the transition mass with smoothing, 
as discussed in Godis et al. (2012). 

Of course, we have here completely ignored the effect 
of feedback, which will play some - yet undefined - role in 
redistributing the cosmic pristine gas falling onto forming 
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galaxies. Another issue would be to estimate for how long 
this entanglement between the large-scale dynamics and the 
kinematic properties of high redshift pervades, given the dis¬ 
ruptions induced by feedback. What will be the effect of 
AGN feedback (Dubois et al. 2013; Prieto et al. 2014) on 
tidally biased secondary infall? Ocvirk et al. (2008) have also 
shown that at lower redshift, the so-called hot mode of ac¬ 
cretion will kick in; how will hot flows wash out/disintegrate 
these ribbons? Given that they locally reflect the large-scale 
geometry, will the gas continue to flow-in along preferred 
directions (as does the dark matter, see e.g. Aubert et al. 
2004), or does the hot phase erase any anisotropy? Will the 
above-mentioned smaller-scale non-linear dynamics eventu¬ 
ally wash out any such trace? 


6.2 Perspectives 

One possibly significant shortcoming of the analysis is the 
proxy involved in using the hessian of the density instead of 
the inertia tensor (though see Section Al). This is critical in 
order to retain a point process for the induced spin, but is 
achieved at the expense of having an adequate estimate for 
the amplitude of the spin, which is unfortunate because from 
the point of view of morphology, the dividing line between 
spirals and ellipticals is likely to be spin amplitude. Let us 
nonetheless assume that e.g. match to simulations or Ansatz 
such as those described in Schafer & Merkel (2012) will yield 
access to reasonable fit to spin amplitude and discuss briefly 
implications to galaxy formation within its cosmic web. 


6.2.1 Epoch of maximal spin advection? 

The inspection of hydrodynamical simulations (e.g. Godis 
et al. 2012, using tracer particles) shows that ribbon-like 
caustics feed the central galaxy along its spin axis from 
both poles. The gas flowing roughly parallel to the spin axis 
of the disc along both directions will typically impact the 
disc’s circum-galactic medium and shock once more (as it 
did when it first reached the wall, and then the filaments, 
forming those above mentioned ribbons), radiating away its 
vertical momentum (see Tillson et al. 2015). These ribbons 
are generated via the same winding/folding process as the 
protogalaxy, and represent the dominant source of secondary 
filamentary infall which feeds the newly formed galaxy with 
gas of well-aligned angular momentum. 

Having computed the most likely spin (direction) as a 
function of position, it is therefore of interest to measure its 
covariant polar flux through a drifting forming galaxy. 

From our knowledge of the spin distribution within the 
neighbourhood of a given saddle we may then compute the 
rate of advected spin within some galactic volume V = SAz] 
it reads 

s = J d^S ‘ V ^ ps = j d^r V • (pv G s) , 

^ S[pw s]+ SAz ^ {pvzSz) , (40) 

where the last equality assumes that the advection is quasi- 
polar, and that the spin is mostly aligned with the filament. 
In equation (40) v is the gradient of the potential. Let us 
identify the curvilinear coordinate, Zup, for which this flux 


is maximal: 


d 

Zup = argmaxT— 

z oz 


{pVzSz) = f 


^ = 0 }. 


(41) 


The coordinate ^up(i^,^ 2 ) characterizes the most active 
regions in the cosmic web for galactic spin up. Focussing 
on the most likely saddle, the argument sketched in Sec¬ 
tion 4.3.2 allows us to assign a redshift-dependent spin-up 
mass, A4up(Z), via equations (37) and (39). There could 
be an observational signature, e.g. in terms of the cosmic 
evolution of the SFR, as maximum spin-up corresponds to 
efficient pristine cold and dense gas accretion, which in turn 
induces consistent and steady star formation. 


6.2.2 Morphological type versus loci on weh? 

The magnitude of the spin of galaxies could be taken as a 
proxy for morphological type. Indeed, Welker et al. (2014, 
2015) have shown in cosmological hydrodynamical simula¬ 
tions that spin direction and galactic sizes where sensitive to 
the anisotropic environment. It is shown in particular that 
the magnitude of the spin of simulated galaxies increases 
steadily and aligns itself preferentially with the nearest fila¬ 
ment when no significant merger occurs, in agreement with 
the first phase of the above described spin-up (see also Pi¬ 
chon et al. 2011). During that phase, the fraction of larger 
spirals should increase. In contrast, following Figure 24, if 
we account for the fact that galactic morphology - the frac¬ 
tion of ellipticals, correlates with dark halo mass, it should 
then increase with distance to saddle. 

In order to tackle such process theoretically, it would 
therefore be worthwhile to revisit Quinn & Binney (1992) in 
the context of this constrained theory of tidal torques and 
quantify how the dynamics of concentric shells are differen¬ 
tially biased by the tides of a saddle point. This would allow 
us to describe the whole timeline of anisotropic secondary 
infall. 


6.2.3 Implication for weak lensing? 

Weak lensing attempts to probe the statistics of the cos¬ 
mic web between background galaxies - which shape is as¬ 
sumed to be uncorrelated - and the observer, while assum¬ 
ing that observed shape statistics reflects the deflection of 
light going through the intervening web. In view of Fig¬ 
ure 12, if we take as a proxy spin alignment for shape 
alignment, we can in principle compute the expectation of 
^(Ar) = (s(r) • s(F)|skl) as a function of Ar = r — rb 
Galling Ss = s — (s(r)|skl), we have ^(Ar) = (s(r)|skl) • 
(s(r')|skl) + (5s(r) • 5s(r')|skl) + 2((5s(r) • s(r')|skl). Let 
us just focus here on the hrst term, (s(r)| ski) • (s(r')| ski). 
Given equation (25), we can compute it and find that it 
will typically be non-zero and vary significantly depending 
on both the magnitude and the orientation of Ar. E.g. if 
Ar is off axis along the filament, but if the pair is close to 
the saddle and |Ar| is small it will be positive (spins will 
align as they are both within coherent region of the saddle’s 
tides), while if |Ar| is somewhat larger it will vanish (spins 
will be perpendicular). Gonversely, if Ar is transverse to the 
filament and |Ar| is small, it will be positive, but if |Ar| is 
of the order of the size of one octant it will typically vanish 
again. The formalism presented in Section 3.1 can clearly be 
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extended (while considering the joint three points statistics) 
to predict exactly all terms involved in ^(Ar) and quantify 
within this framework the effect of intrinsic alignments on 
the spin-spin two point correlation. This is will be to topic 
of future work (Codis et al, in prep.). 


ACKNOWLEDGMENTS 

This work is partially supported by the Spin(e) grants ANR- 
13-BS05-0005 (http://cosmicorigin.org) of the French 
Agence Nationale de la Recherche and by the ILP LABEX 
(under reference ANR-lO-LABX-63 and ANR-ll-IDEX- 
0004-02). CP thanks D. Lynden-bell for suggesting to tackle 
this problem and Churchill college for hospitality while 
this work was completed. Many thanks to J. Devriendt, 
A. Slyz, J. Binney, Y. Dubois, V. Desjacques, C. Laigle 
and S. Prunet for discussions about tidal torque theory, 
and to our collaborators of the Horizon project (http:// 
projet-horizon.fr) for helping us produce the Horizon-47r 
simulation, E. Bouchet for allowing us to use the magiqueS 
supercomputer during commissioning, and to S. Rouberol 
for making it possible, and running the horizon cluster for 
us. SC and CP thanks Lena for her hospitality when this 
work was initiated, and Eric for his help with some figures. 


References 

Aragon-Calvo M. A., 2013, preprint (arXiv: 1303.1590) 
Aragon-Calvo M. A., van de Weygaert R., Jones B. J. T., van der 
Hulst J. M., 2007, ApJ Let, 655, L5 
Aubert D., Pichon C., Colombi S., 2004, MNRAS, 352, 376 
Bailin J., Steinmetz M., 2005, ApJ, 627 

Bardeen J. M., Bond J. R., Kaiser N., Szalay A. S., 1986, ApJ, 
304, 15 

Bernstein G. M., Norberg R, 2002, AJ, 124, 733 
Bond J. R., Cole S., Efstathiou G., Kaiser N., 1991, ApJ, 379, 
440 

Bond J. R., Kofman L., Pogosyan D., 1996, Nature, 380, 603 
Brown M. L., Taylor A. N., Hambly N. C., Dye S., 2002, MNRAS, 
333, 501 

Catelan P., Theuns T., 1996, MNRAS, 282 

Catelan P., Kamionkowski M., Blandford R. D., 2001, MNRAS, 
320, L7 

Codis S., Pichon C., Devriendt J., Slyz A., Pogosyan D., Dubois 
Y., Sousbie T., 2012, MNRAS, 427, 3320 
Codis S., et al., 2015, MNRAS, 448, 3391 

Crittenden R. G., Natarajan P., Pen U.-L., Theuns T., 2001, ApJ, 
559, 552 

Croft R. A. C., Metzler C. A., 2000, ApJ, 545, 561 
Danovich M., Dekel A., Hahn O., Teyssier R., 2012, MNRAS, 
422, 1732 

Doroshkevich A. G., 1970, Astrophysics, 6, 320 
Dubois Y., Pichon C., Haehnelt M., Kimm T., Slyz A., Devriendt 
J., Pogosyan D., 2012, MNRAS, 423, 3616 
Dubois Y., Pichon C., Devriendt J., Silk J., Haehnelt M., Kimm 
T., Slyz A., 2013, MNRAS, 428, 2885 
Dubois Y., et al., 2014, MNRAS, 444, 1453 
Forero-Romero J. E., Contreras S., Padilla N., 2014 
Gunn J. E., Gott HI J. R., 1972, ApJ, 176, 1 
Hahn O., Porciani C., Carollo C. M., Dekel A., 2007, MNRAS, 
375, 489 

Heavens A., Refregier A., Hey mans C., 2000, MNRAS, 319, 649 
Hey mans C., Brown M., Heavens A., Meisenheimer K., Taylor A., 
Wolf C., 2004, MNRAS, 347, 895 


Hirata C. M., Seljak U., 2004, Rhys. Rev. D, 70, 063526 
Hirata C. M., Mandelbaum R., Seljak U., et al. 2004, MNRAS, 
353, 529 

Hirata C. M., Mandelbaum R., Ishak M., Seljak U., Nichol R., 
Pimbblet K. A., Ross N. P., Wake D., 2007, MNRAS, 381, 
1197 

Hoyle F., 1949, Problems of Cosmical Aerodynamics, Central Air 
Documents, Office, Dayton, OH. Central Air Documents Of¬ 
fice, Dayton, OH 

Huchra J. R, Geller M. J., 1982, ApJ, 257, 423 
Joachimi B., Mandelbaum R., Abdalla F. B., Bridle S. L., 2011, 
A&A, 527, A26 

Joachimi B., Semboloni E., Bett P. E., Hartlap J., Hilbert S., 
Hoekstra H., Schneider P., Schrabback T., 2013a, MNRAS, 
431, 477 

Joachimi B., Semboloni E., Hilbert S., Bett P. E., Hartlap J., 
Hoekstra H., Schneider R, 2013b, MNRAS, 436, 819 
Klypin A., Shandarin S. F., 1993, ApJ, 413, 48 
Laigle C., et al., 2015, MNRAS, 446, 2744 
Lee J., Pen U., 2000, ApJ, 532, L5 
Lee J., Pen U.-L., 2002, ApJ Let, 567, Llll 
Libeskind N. 1., Hoffman Y., Forero-Romero J., Gottlober S., 
Knebe A., Steinmetz M., Klypin A., 2013, MNRAS, 428, 2489 
Mandelbaum R., Hirata C. M., Ishak M., Seljak U., Brinkmann 
J., 2006, MNRAS, 367, 611 

Mandelbaum R., Blake C., Bridle et al. 2011, MNRAS, 410, 844 
Musso M., Sheth R. K., 2012, MNRAS, 423, L102 
Ocvirk P., Pichon C., Teyssier R., 2008, MNRAS, 390, 1326 
Paranjape A., Lam T. Y., Sheth R. K., 2012, MNRAS, 420, 1429 
Paz D. J., Stasyszyn F., Padilla N. D., 2008, MNRAS, 389, 1127P 
Peacock J. A., Heavens A. F., 1990, MNRAS, 243, 133 
Peebles P. J. E., 1969, ApJ, 155, 393 
Pichon C., Bernardeau F., 1999, A&A, 343, 663 
Pichon C., Pogosyan D., Kimm T., Slyz A., Devriendt J., Dubois 
Y., 2011, MNRAS, pp 1739-+ 

Pogosyan D., Bond J. R., Kofman L., Wadsley J., 1998, in 

S. Colombi, Y. Mellier, & B. Raban ed.. Wide Field Surveys 
in Cosmology, p. 61 (arXiv: astro-ph/9810072) 

Pogosyan D., Pichon C., Gay C., Prunet S., Cardoso J. F., Sousbie 

T. , Colombi S., 2009, MNRAS, 396, 635 

Porciani C., Dekel A., Hoffman Y., 2002, MNRAS, 332, 325 
Prieto J., Jimenez R., Haiman Z., Gonzalez R. E., 2014, preprint 
(arXiv:1410.4061) 

Prunet S., Pichon C., Aubert D., Pogosyan D., Teyssier R., Got- 
tloeber S., 2008, ApJ Sup., 178, 179 
Quinn T., Binney J., 1992, MNRAS, 255, 729 
Schaefer B. M., 2009, International Journal of Modern Physics D, 
18, 173 

Schafer B. M., Merkel P. M., 2012, MNRAS, 421, 2751 
Schneider M. D., Bridle S., 2010, MNRAS, 402, 2127 
Schneider M. D., Frenk C. S., Cole S., 2012, JCAP , 5, 30 
Sousbie T., Pichon C., Colombi S., Pogosyan D., 2008, MNRAS, 
383, 1655 

Sousbie T., Colombi S., Pichon C., 2009, MNRAS, 393, 457 
Spergel D. N., et al., 2003, ApJ Sup., 148, 175 
Stewart K. R., Brooks A. M., Bullock J. S., Mailer A. H., Diemand 
J., Wadsley J., Moustakas L. A., 2013, ApJ, 769, 74 
Tenneti A., Singh S., Mandelbaum R., Matteo T. D., Feng Y., 
Khandai N., 2015, MNRAS, 448, 3522 
Teyssier R., 2002, A&A, 385, 337 
Teyssier R., et al., 2009, A&A, 497, 335 

Tillson H., Devriendt J., Slyz A., Miller L., Pichon C., 2015, MN¬ 
RAS, 449, 4363 

Welker C., Devriendt J., Dubois Y., Pichon C., Peirani S., 2014, 
MNRAS, 445, L46 

Welker C., Dubois Y., Devriendt J., Pichon C., 2015, preprint 
(arXiv:1510.4061) 

White S. D. M., 1984, ApJ, 286, 38 


MNRAS 000 , 000-000 (0000) 


22 S. Codis, C. Pichon and D. Pogosyan 


White S. D. M., Tully R. B., Davis M., 1988, ApJ Let., 333, L45 
Zel’dovich Y. B., 1970, A&A, 5, 84 

Zhang Y., Yang X., Faltenbacher A., Springel V., Lin W., Wang 
H., 2009, ApJ, 706, 747 


APPENDIX A: A MULTI-SCALE THEORY 

The proxy we take for the spin direction 

Si — ^ ^ ^ijkHjiTiJ^ , (Al) 

is, as mentioned in the main text, a (quadratic) approxima¬ 
tion. First, because Equation (1) is only valid in the linear 
regime (Porciani et al. 2002), but possibly more importantly 
because we take the Hessian as a proxy for the inertia ten¬ 
sor. In practice, recall that this approximation seems nev¬ 
ertheless to capture the essence of the processes at work in 
aligning spins with the large-scale structure given its abil¬ 
ity to explain observed alignments through the comparison 
with simulations presented in Section 5.2 together with the 
measured mass transition. This suggests experimentally that 
it is indeed reasonable. Notwithstanding, while Hij and lij 
locally share the same eigenframe, their amplitudes are dif¬ 
ferent, leading to a different weighting of field configurations 
when computing ensemble averages such as in equation (14). 
It is therefore important to investigate this possible short¬ 
coming further in this Appendix. 


Al More realistic spin proxies and peak 


For this purpose, one can in principle i) impose an additional 
peak constraint at the location where the spin is computed 
in order to impose that a proto-halo will form there, and ii) 
use more realistic spin proxies. 

The peak constraint will typically be at a smaller scale 
than the filament’s constraint, which requires building a two- 
scale theory and therefore increases significantly the com¬ 
plexity of the formalism. An additional difficulty with ii) is 
that standard local proxies for the inertia tensor are highly 
non-linear and therefore require high-dimensional numerical 
integrations that are fairly difficult to implement in practice. 
For instance, considering the proxy that Schafer & Merkel 
(2012) use to locally approximates the inertia tensor, we 
have 


lij — 


M 

y 


aIj-aI 0 0 

0 0 
0 0 A^ T A^ 


(A2) 


(in the frame of the Hessian) where the mass is M = 
AjSTrAxAyAzPoaQ and the semi-axes of the ellipsoid, A^, are 
function of the eigenvalues of the Hessian (negative for a 
peak). 


,A3) 

The traceless part of lij that is relevant for torques is then 
proportional to the traceless part of the inverse Hessian 

= ‘^va2MHj . (A4) 

This introduces singular factors like 1 / Vdet H in the expec¬ 
tation for (s). Such factors make the numerical evaluation of 



Log RhIRs 


Figure Al. Maximum of spin as a function of the halo’s scale 
Rh- 


equation (24) more challenging as discussed in Section 2. We 
therefore postpone their evaluation in three dimensions to 
future work. Let us briefly investigate their implementation 
in two dimensions. In Section A2 we introduce a multi-scale 
description while in Section A3 we take into account the 
proxy given by equation (A2) and we add an explicit peak 
condition. 


A2 A 2D multiscale analysis 


The evaluation of r*, see equation (20), requires to look for 
the scale that maximizes the spin amplitude. An improve¬ 
ment to the main text’s approach is to investigate this issue 
in a two-scale theory where the scale of the central peak is 
larger than the scale of the halo. In this section, we propose 
to describe how to implement this multi-scale theory and 
we show that this more accurate multi-scale estimate for r* 
yields, following the main text, a very similar value for the 
critical mass. 

More specifically, let us smooth the fields in yy (the 
location of the spin) on Rh^ an additional parameter, char¬ 
acterizing the halo’s size and the fields in yx on a different 
scale Rs > Rh in order to impose a large-scale filament (a 
peak in 2D). If the one-point covariance matrix Co, does not 
change for a power law power spectrum, the two-point co- 
variance matrix does. In particular, all correlation func¬ 
tions ^ are now function of r, Rs and Rh- For instance for a 
power-law density power spectrum, P{k) oc kA, 


M = 


Rl + Rl ^ 

2RsRh 


T?-lin + 2) ^ 


Rl + Rl 


^2 


where Tl — iFi (n/2 + z; j; —r^/2(Rs + R^)). The mean 
spin in yy is then given by the expectation of Sz — 
'Yhi j k ^i 3 ^ 4 >ikXjk given a peak on scale Rs in yx as was com¬ 
puted in the main text. Compared to the main text, the only 
difference here is that we now also take into account the two- 
scale process through the two smoothing scales, Rs and Rh- 
The maximum spin magnitude as a function of the scale 
Rs/Rh is then computed and displayed in Figure Al. It ap¬ 
pears that the spin magnitude is non-monotonic, peaking at 
Rh = O.SRs which is very close to the value of r* ~ 0.7Rs 
(when top hat smoothing are taken for both lengths). 

Note that as we are computing here only the component 
of the spin along the filament, our magnitude plot does not 
contain the mass prefactor, and our proxy for the moment 
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of inertia, I H, reflects only its orientation, but not its 
magnitude. We are able to argue, however, that Figure A1 
displays a maximum spin alignment with the filament’s di¬ 
rection for some critical value of Rh - as it shows first an 
increase and then a fall in the z-component of the spin that 
is due to Hessian-tidal shear alignment. 


A3 A 2D multiscale analysis with peak constraint 


Adding a peak constraint at the location of the spin and tak¬ 
ing into account the proxy given by equation (A 2 ) is more 
difficult as it requires a numerical integration to account for 
the sign constraints on the eigenvalues. Notwithstanding this 
shortcoming, we will show now that in two dimensions, the 
adjunction of a peak constraint preserves both the qualita¬ 
tive picture (same geometry with four quadrants of opposite 
spin direction) as well as the typical scale for r*. 

Figure A2 indeed shows the numerical integration of 
the mean spin when i) the inertia tensor is approximated by 
equation (A4) -where the mass is fixed by the smoothing 
length Rh-^ ii) the scale of the central peak Rs is different 
from the scale of the halo Rh (here we take Rh/Rs = 1 / 10 )? 
iii) there is a peak constraint at the location where the spin 
is computed with height u = 5/2, negative eigenvalues and 
zero gradient. In short, the mean spin is now computed as 


(s|pk,pk) 


(sdet[H] 0(-Ai)5D(a: - !^)(5D(a:i)|pk) 
(det[H]0(-Ai)<5D(a:-O<5D(a:i)|pk) ’ ^ ’ 


where s, is defined as s, = ^ ^ tijkHJj^Tik while the ex- 

pectations ('Ipk) are defined as conditional expectation to 
a central peak of geometry = 1, Ai = —1, A 2 = —2. The 
mean spin map is then obtained by numerical integration. 
Figure A 2 clearly shows that the four quadrants of oppo¬ 
site spin direction, as well as the size of these quadrants are 
preserved. This test strongly suggests that in two dimen¬ 
sions, improvements beyond the I H approximation do 
not change the global picture described in the main text. 

Further developments, beyond the scope of this paper, 
could be to carry out the same analysis in three dimensions, 
also adding a peak constraint at the location where the spin 
is computed in order to impose the existence of a proto-halo 
and use equation (A 2 ) to define its inertia tensor. While the 
two-scale analysis is straightforward enough to implement, 
the adjunction of a peak constraint in three dimensions is 
much more tricky and requires in particular the computation 
of high-dimension numerical integrals (the results will not be 
analytic anymore) that are left for future investigations. 


APPENDIX B: DUAL VOID THEORY 

The theory presented in Section 4.1 is algebraic. Effectively 
no assumption has been made about the signs of the eigen¬ 
values of the saddle we are considering. It is therefore also 
perfectly valid in the neighbourhood of a wall-type saddle 
in order to describe the spin alignments of dark haloes in 
that vicinity. At a qualitative level. Figure 3 applies up to a 
sign: voids and wall saddles repel. It follows that the spins 
should rotate around the wall-saddle to void axis and be¬ 
come parallel near the wall with a point symmetric change 
of polarity. This is indeed what equation (25) predicts and 
is shown on Figure BI. 



x/Rs 


Figure A2. Mean spin computed via numerical integrations 
when i) a central peak with geometry given hy u = I, Ai = 
— I, A 2 = —2 is set, ii) the inertia tensor is approximated by 
equation (A2), iii) the scale of the central peak Rs is different 
from the scale of the halo Rh (here we take Rh/Rs = I/IO)? 
iv) there is a peak constraint at the location where the spin is 
computed with height u = 5/2, negative eigenvalues and zero 
gradient. The density power spectrum is a power-law here with 
spectral index n = 1/2. 


The statistical significance of these alignments is likely 
to be reduced as there are much fewer galaxies in voids and 
near wall saddles. 


APPENDIX C: TECHNICAL COMPLEMENTS 

Cl Codes for density and spin in 2/3D 

The expression for the 2 and 3D spin statistics (mean 
and variance) for scale invariant power spectra are avail¬ 
able both as a mathematica package (http://www.iap. 
fr/users/pichon/spin/code/ATTT .m), and a mathemat¬ 
ica notebook (http://www.iap.fr/users/pichon/spin/ 
code/ATTT-package .nb). The following functions are pro¬ 
vided: (52D, spin2D, var2D, (53D, spinSD, which corre¬ 
spond resp. to the 2D density, spin, its variance, and in 3D 
the density and the spin for scale invariant power spectra of 
index n as a function of position r, 0, ( 0 ) and the geometry 
of the peak (resp. saddle) Ai, A2, (A3). Compiled versions 
are also provided. 


C2 Correlation functions for power-law spectra 

The 2D correlation functions defined in Equation (5) and 
(9) can be analytically obtained for density power-law power 
spectrum Pk{k) oc k'^ for which the scale parameter is 7 = 
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X / Rs X / Rs 

Figure Bl. Top left-hand panel: Mean density in the plane of the wall Oxy and centered on a wall-type filament with geometry u = 0.5, 
Ai = 0.8 along the x-axis, A 2 = 0.6 along the y-axis and A 3 = —0.5 along the z-axis for a power-law density power spectrum with spectral 
index n = — 2 . Contours are displayed from S = 0.6 to 1.2 as labeled. The filaments are clearly seen around y = ± 2 i^s. Top right-hand 
panel: Mean spin colour-coded by its projection along the normal to the wall. The spins are aligned with the normal in the plane of the 
wall and perpendicular to it when going outside the plane of the wall. Bottom panels: mean spin vectors at 2 : = 0.5i^s (left) and 2 : = IRg 
(right). Contours represent the orientation of the spin with regards to the normal to the wall from -1 (anti-aligned, red) to +1 (aligned, 
blue) through 0 (perpendicular). 


y^(?Tr+~2)7(?7H-^: 

<^AA/ \ /^ (^ + 4 ) 2 ^3 7-2^ 

= 7 I IQ ^ 

^x^{r) = 

e^/W=-7(n + 2)r2j-|, 
^{n+Ayxi 
y^{r)=-^{n + Q)pTl 


iy{r)=^-{AXl-ZFl), (Cl) 

where r is in units of the smoothing length and = 
iFi (n/2 + z; —r^/4), with iFi the Hypergeometric func¬ 
tions of the first kind. Some of those correlation functions 
are plotted in Figure 5. 

The 3D correlation functions can similarly be obtained 
for power-law power spectrum Pk{k) oc k'^{soine of those 
correlations are plotted in Figure 10) for which the scale 
parameter is 7 = y^(n + 3)/(n + 5) and we define Gj =1 
Fi (n + i/2;j72;-rV4): 
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st^w= 




&A(r) = 


-32(n - 1)03 _ ((2 - 4n)r2 + - 32) ^li + (n - 2) (2nr2 - + 32) 

2(n^ — 1) ’ 

r (^) (32(n + l)gg + (-2(2n + 3)r^ + - 32) + n (-2(n + 2)r2 + - 32) Qf) 

-32(n + 3)6| + (2(2n + 7)r^ -A + 32) 6^ + (n + 2) (2(n + A)r^ -A + 32) gf 

2(n + 3)(n + 5)r2 ’ 

(4n (r2 + 3) - (r^ + 8) - 40) + 16(n - l)gf + (n - 2) (2(n - 3)r'^ - - 28) gJii 

2 {v? — l)r^ ’ 

r (^) (-16(n + 1)6| + (-4n {r^ + 3) + + 16) g} + n (-2(n - l)r2 + + 28) gf) 




16(n + 3)^1 + (4(n + 2)r^ + 12n — + 8 ) ^3 + (n + 2) (2(n + l)r^ — — 28) ^3 

2(n + 3)(n + 5)r2 

((n - 2)r2 p + 10) - 48) + (-2(n - 6)r2 + + 48) 


(n^ — 1) 

r (^) ((nr2 p + 10) - 48) Gf + (-2(n - 4)r‘^ + + 48) Gl) 

((n + 2)r^ (r^ + lO) — 48) G^ + (—2(n — 2)r^ + + 48) G\ 

{n + 3)(n + 5)r^ 



Figure Cl. Two-point correlation functions as a function of the 
separation r in units of the smoothing length Rg = 5Mpc//i for a 
ACDM power spectrum. 


C3 Correlation functions for LCDM spectra 

The same ^ correlation functions can also be computed for 
a ACDM power spectrum using Bardeen et al. (1986) and 
equation (9). The corresponding functions are shown on Fig¬ 
ure Cl for a Gaussian smoothing length of Rg = 5Mpc//i 
and a WMAP-7 cosmology. Note that those correlation func¬ 
tions are quite similar to n = —2 power-law power spec¬ 
trum (see Figure 10). Given these correlations, it would be 
straightforward to compute the corresponding spin. 
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